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Abstract 

Quasielastic electromagnetic and parity- violating electron scattering response 
functions of relativistic nuclear matter are reviewed. The roles played by the 
Hartree-Fock field and by nuclear correlations in the Random Phase Approx- 
imation (treated within the continued fraction scheme) are illustrated. The 
parity-violating responses of nuclei to polarized electrons are also revisited, 
stressing in particular the crucial role played by the pion in the nuclear dy- 
namics. Finally, some issues surrounding scaling and sum rules are addressed. 
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I. INTRODUCTION 



Traditionally, studies of excitations of the nucleus via the interactions between leptons 
and the nucleus have been centered mostly on the longitudinal (Rl) and transverse (-Rt) 
nuclear response functions explored in unpolarized, inelastic, inclusive (e, e') electron scat- 
tering. The present paper will also be largely restricted to the investigation of these quan- 
tities, although the set of experimentally accessible responses is much larger, including also 
semi-inclusive responses and responses that arise when initial or final hadronic degrees of 
freedom are active. In the former, particles are detected in coincidence with the scat- 
tered electron ((e, e'p), (e,e'n), (e,e'pn), etc.), whereas the latter includes reactions such as 
A(e,e'), A(e,e'p ), etc. Of course, in all cases the incident electron may also be polarized. 
A concerted effort is presently being placed on experimental studies of this extended set of 
responses where special sensitivities to normally hidden aspects of nuclear structure are ex- 
pected to exist. Such studies are very challenging and thus only relatively recently have they 
been made feasible by advances in accelerators, in developments of polarized electron beams 
and polarized nuclear targets, and in the construction of the required hadron polarimetry. 
Issues still surround the inclusive unpolarized responses themselves, however, and since a 
successful level of understanding of the underlying nuclear dynamics would be incomplete 
without a coherent picture of the entire set of responses, unpolarized inclusive and semi- 
inclusive/polarized, the former still deserve continued study and so provide the focus for the 
present article. 

Beyond the electromagnetic (EM: parity-conserving, vector) responses i?L and -Rt our 
study will also include their parity-violating analogs R^ v and R^ as well as the nuclear 
parity- violating axial response R^t". Here AV indicates that the axial leptonic and vector 
hadronic currents enter; VA indicates the converse. This larger set of inclusive responses may 
be explored via the inclusive scattering of longitudinally polarized electrons from unpolarized 
nuclei, since the electron helicity asymmetry is parity-violating. The three new responses 
all arise from interferences between the weak neutral current (WNC) and EM current. In 
the cases of Rf^ and i?^ V it is the vector part of the WNC that enters and this is believed 
to be closely related to the EM current (in the absence of strangeness content in the nucleus 
these two responses are tied to R^ and -Rt); in the case of R^t" it is the axial part of the 
WNC that enters, namely an interesting new inclusive nuclear response function. 

The vector responses (the four labeled either L or T) are, of course, interesting in their 
own right, since disentangling them through measurements of both parity-conserving and - 
violating inclusive electron scattering would permit the isolation of the isoscalar and isovector 
contributions they contain, as will be discussed in detail later. Accomplishing this separation 
would represent a significant step forward in our understanding of nuclear structure: indeed, 
for a long time researchers have sought possible ways of "measuring" how nuclear correlations 
work in the isoscalar and isovector channels. 

A further point worth noting is that analogous responses also play a role in the scattering 
of hadrons from nuclei. However, in order to interpret the experimental data properly, in 
addition to the response functions there one also needs a reliable description of the reaction 
mechanism. In fact, unlike either real or virtual photons, hadrons are mostly absorbed or 
scattered at the surface of the nucleus. Moreover, the hadrons disrupt the nucleus to a 
much larger extent than do photons, and therefore the interpretation of reactions induced 
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by photons is generally felt to be under better control than those induced by hadrons. 
Moderating this statement to some degree is the fact that electron scattering is still somewhat 
flawed by a not yet fully satisfactory understanding of dispersive effects and of the distortion 
of the electron waves moving in the EM potential of the nucleus. Accurate knowledge of the 
nuclear response functions gained with electron scattering provides a way to test the reaction 
mechanisms of the models employed in hadron scattering. However, the electroweak studies 
do not provide all of the information we seek and in this regard it is worth noting that in 
some cases reactions induced by hadrons give access to nuclear responses that are not easily 
extracted in electron scattering, the best example in this connection being offered by the 
long sought after spin-longitudinal isovector response. 

Turning now to the problem of modeling the nuclear response functions, in the present 
work we confine our attention to the quasifree region, which is well suited for a microscopic 
treatment in terms of nucleons and mesons, specifically in terms of the standard field the- 
oretical techniques that we employ. Although the A peak may also be treated in the same 
framework, it will not be dealt with here to curb the length of this article. Since we are 
concerned with kinematical regions where relativistic effects are relevant, a good starting 
point for the development of our approach is given by the Relativistic Fermi Gas (RFG), a 
covariant model in the sense that its ingredients are the fully relativistic nucleon propagators 
and EM/WNC vertices. Of course the RFG misses surface and finite-size effects. First of 
all, these are of secondary importance in obtaining a general understanding of the scattering 
of electrons in the quasielastic and A peak domains. Secondly, they can be accounted for 
within the semiclassical approach, which exploits the advantages offered by the translational 
invariance of the RFG and yet is able to incorporate some of the physics of a finite system. 

As discussed in detail later, the perturbative approach we follow requires the setting up 
of a nuclear mean field (Hartree-Fock, HF) and the treatment of the residual interaction 
effects in the fully antisymmetrized Random Phase Approximation (RPA). In addition one 
needs as preliminary input the nucleon-nucleon force: since we would like to view the nucleus 
as an interacting system of baryons and mesons, a natural choice in this connection is given 
by a meson-exchange interaction such as the Bonn potential, which can be cast in the 
framework of an effective field theory. Another preliminary problem relates to the short- 
range nuclear correlations induced by the violent repulsion present in the nucleon-nucleon 
force at small distances. A technique for their treatment is indeed available, namely the 
summation of the Brueckner ladder diagrams; however, it is not yet possible to employ this 
technique covariantly especially at high density, where on the one hand ladder diagrams are 
increasingly important and on the other the role of relativity cannot be ignored. In lieu of 
this, in a few of the results discussed in the following section we shall indicate what insight 
can be gained by employing a parameterization of a non-relativistic G-matrix based on the 
Bonn potential. 

II. QUASIELASTIC RESPONSE FUNCTIONS FOR INCLUSIVE ELECTRON 

SCATTERING 

As mentioned in the Introduction, in past years quasielastic electron scattering from 
nuclei has been the subject of intense experimental [1-3] and theoretical (see, e. g., Refs. [4]- 
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[18]) investigations. The first aim of the theoretical studies is to test the available nuclear 
models; once the nuclear physics issues are well understood, one might then hope to gain 
insight into other aspects of the problem, for instance into the form factors of the nucleon, 
which can be extracted from the data with an accuracy that is strictly connected to our 
ability to handle the nuclear physics. 

In principle, the quasifree regime is thought to be the obvious place to focus on, as 
one hopes that there the physical quantities of interest may be computed in a reliable way, 
while also in this case in practice one has to cope with significant problems. Many diverse 
techniques have been employed in the literature. Each of them has its own relative merits 
and deficiencies and clearly it would be highly desirable to be able to reach some degree of 
convergence in their outcomes. 

In the following [19], we shall be concerned with Green's function techniques as intro- 
duced, e. g., in Ref. [20]. This method can be, and has been, applied both to finite nuclei 
and nuclear matter. Here, we shall focus on nuclear matter, having in mind applications 
to electron scattering (that is, without the complications introduced by the reaction mech- 
anism of hadronic probes) in a range from a few hundreds to about 1 GeV/c of transferred 
momentum where the quasielastic peak is far from low-energy resonances and not too much 
affected by finite-size effects. The use of nuclear matter reduces the computational load, thus 
allowing a more straightforward implementation of more sophisticated theoretical schemes 
than would otherwise be feasible, and this makes it easier to develop and test approximation 
methods that might subsequently also be utilized for calculations in finite nuclei. 

Let us now briefly summarize the theoretical framework that we shall discuss in detail 
in the following subsections. 

A first issue one has to confront in setting up the formalism concerns the treatment of 
relativistic effects. Kinematical effects, while obviously rather important, can be included in 
a straightforward way. The treatment of dynamical effects is more delicate. Two main paths 
have been followed in the literature, either using field theoretical methods (as done, e. g., in 
the Walecka model and its derivations [21]) or using potential techniques (i. e., employing 
phenomenological potentials truncated at some order in the non- relativistic expansion). Here 
we shall take the second path, but to limit the amount of material to be covered, we shall 
discuss only non-relativistic potentials. 

The extensions necessary to include higher-order relativistic terms are discussed in 
Ref. [14], where the nuclear response functions have been calculated using techniques similar 
to the ones explained below, using as an input the relativistic Bonn potential [22] expanded 
in powers of P/m N and q/m N up to second order — P and q being the average of the incom- 
ing and outgoing nucleon momenta and the exchanged momentum, respectively. As shown 
in [14], the effect of these dynamical relativistic corrections is significant; indeed, the valid- 
ity of that expansion at high momenta and the inclusion in that framework of short-range 
nucleon-nucleon correlations has yet to be explored (see, however, Refs. [14,23,24]). 

Next, one should choose the phenomenological input potential and, in connection with 
this choice, attempt to cope with the problem of dealing with short-range correlations. All 
of the formulae given in the following sections are based on a generic one-boson-exchange 
potential. They can thus be used both with a bare phenomenological interaction — such as 
one of the Bonn potential variants — or with a one-boson-exchange parameterization of the 
G-matrix generated from some potential. The use of an effective interaction derived from a 
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G-matrix is a common way of including short-range correlations. However, apart from the 
relativistic issue, one should be aware of possible problems due to the use of a local potential 
to fit non-local matrix elements. At least in a few cases discussed in the literature this does 
not appear to be a reason for concern [25,26]. On the other hand, possible effects arising 
only in the quasielastic regime remain completely unexplored. Indeed, G-matrices employed 
in quasielastic calculations are usually generated using bound-state boundary conditions, 
which make them real and practically energy- independent, while in general they are both 
complex and energy-dependent. 

Once we have fixed the effective interaction, we can proceed to consider a hierarchy of 
approximation schemes. 

The lowest-order approximation is, of course, given by the free Fermi gas. Then, one 
may include mean-field correlations at the HF level (or Brueckner-Hartree-Fock (BHF) if 
short-range correlations are accounted for). In nuclear matter a HF calculation can be done 
exactly without too much effort. Later we show how a quite accurate analytic approximation 
can be derived, and how this is needed to combine the HF and RPA schemes. The latter is 
the last resummation technique we shall discuss. It should be noticed that even in nuclear 
matter the calculation of the antisymmetrized RPA response functions is not trivial. Indeed, 
most calculations, labeled "RPA" in the literature, are actually performed in the so-called 
"ring approximation", where only the direct contributions are kept. For this case, in nuclear 
matter one gets a simple algebraic equation for the response. Here, we use the continued 
fraction (CF) technique to provide a semi-analytical estimate of the full RPA response 
(see Refs. [6] and [8] for alternative methods). Calculations with this method have been 
performed both in finite nuclei [4,5] and in nuclear matter [27,28,14,23], always truncating the 
CF expansion at first order because of the difficulty of the numerical calculations involved. 
We have pushed the analytical calculation far enough to yield not only a fast and accurate 
estimate of the first-order CF contribution, but also of the second-order one. Since the rate 
of convergence of the CF expansion cannot be assessed on the basis of general theorems, 
this is the only way of getting a quantitative grip on the quality of the approximation. 
As mentioned before, HF (and kinematical relativistic) effects can then be incorporated in 
the RPA calculation, yielding as the final approximation scheme a HF-RPA (or BHF-RPA) 
response function. 

Of course, several many-body contributions have been left out in our analysis. However 
the classes of many-body diagrams discussed here already allow one to capture the main 
features of the quasielastic response and, since semi-analytical methods have been developed 
for their computation, our formalism constitutes a valid starting point for the study of other 
many-body effects. 

A. Response functions 

We consider an infinite system of interacting nucleons at some density fixed by the Fermi 
momentum For the kinetic energies of the nucleons we can choose either relativistic or 
non-relativistic expressions, whereas we assume that the interactions take place through a 
non-relativistic potential. For the latter the following expression in momentum space is 
assumed 



5 



V(k) = V (k) + V T (k) Tl • r 2 + V a (k)<r x ■ <r 2 + V aT (k)<Ti • <x 2 n • r 2 

+V t (k)S 12 (k) + V tT (k)S 12 (k)r 1 • r 2 , (2.1) 



where Si 2 is the standard tensor operator and V a {k) represents the momentum space poten- 
tial in channel a. Here V a (k) has the general form of a static one-boson-exchange potential 
so that in each spin-isospin channel, namely (0, r, a, err, t, tr), it is represented as a sum of 
contributions from different mesons, V a = J2iV^\ In the central channels (0, r, er, err) the 
contribution from any meson can be expressed as the combination of a short-range ("<5") 
piece and a longer range ("momentum-dependent") piece 1 : 

v ^^iW) t <2 - 2a) 

«W=A^(|^)', « = 0,1,2, (2.2b) 
whereas in the tensor channels (t, tr) one has 

^W^a^(f^)', « = 0,1,2. (2.2c) 

In Eqs. (2.2), gf \ g$ D and are the (dimensional) coupling constants of the i-th meson, 
rrii is its mass and Aj the cut-off; more generally, potentials without form factors or with 
monopole or dipole form factors are allowed. 

Our starting point [29-31] is the Galitskii-Migdal integral equation for the particle-hole 
(ph) four-point Green's function 2 , 

G% a& (K + Q,K;P + Q,P) = -G ai (P + Q) G 5p {P) (2rt) 4 5(K - P) 
+iG aX (K + Q)G x ,p(K) J 0- 4 TZ,^(K + Q,K;T + Q,T) G^,^ 5 (T + Q,T; P + Q, P), 

(2.3) 

diagrammatically illustrated in Fig. 1. In Eq. (2.3), G represents the exact one-body Green's 
function, whereas T 13 is the irreducible vertex function in the ph channel. 
Given G ph one can then define the polarization propagator 

II a /3 i7< 5(Q) = n a/ 3 j7( $(g,u;) 

= * / ^0r^ {K +Q> K -> P +Q>n (2.4) 



The nomenclature stems from the fact that, in the absence of form factors, V$ is a constant and 
is represented by a Dirac (5-function in coordinate space, whereas Vmd is the momentum-dependent 
piece. 

2 Capital letters refer to four-vectors and lower-case letters to three-vectors; the Greek letters 
a, (3, ... refer to a set of spin-isospin quantum numbers. 
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FIG. 1. Diagrammatic representation of the Galitskii-Migdal integral equation for the ph 
Green's function, G ph ; T 13 is the irreducible vertex function in the ph channel; the heavy lines 
represent the exact one-body Green's functions. 



whose diagrammatic representation is displayed in Fig. 2. Note that for 11(5,0;) one cannot 
in general write down an integral (or algebraic) equation. 

In the case of electron scattering, one can define charge — or longitudinal — and magnetic 
- or transverse — polarization propagators. These, in the non-relativistic regime, read 

nl(q,u)=tT[dllL(q,u)dl] (2.5a) 
l£(q,u) = J2^K(q,"), 4(g,o0 = tr[6^n(g,^)6y (2.5b) 

ij 

Kj = (Sij - qiqj/q 2 )/2, 



where, for brevity, the dependence upon the spin-isospin indices has been represented in 
matrix form, introducing hats to indicate matrices. In Eqs. (2.5), / labels the isospin 
channel and the longitudinal and transverse vertex operators are given by: 

Ot° = 1/2 / Oft = tn/2 

Ol =1 = n/2 \ 6^ = an/2. ^ 

The inelastic inclusive scattering cross section where the momentum q and energy 00 are 
transferred to the nucleus is a linear combination of the imaginary parts of H^(q, u). It is 
then customary to define longitudinal and transverse response functions according to 

R L , T (q, = R$(q, + R$(q, u), (2.7) 

which are related to IIl,t by 



7 




RLt(v,u) = -—fl%'(q,u})hnn{ T (q,u) 



7T 

3nA f (J) 2 



2A; 



(g,c)Imn[ iT (g,cu), (2.8) 



(T\ 2 

where V is the volume, A the mass number and the /f^ embody the squared EM form 
factors of the nucleon. The latter are briefly discussed in Appendix A. 



B. Non-relativistic vs relativistic kinematics 

The response functions introduced above have been defined as functions of the momen- 
tum transfer q and energy transfer uj. Actually, it is possible — and convenient — to define a 
scaling variable ip that is a function of q and uj and may be used in place of uj. This variable 
is such that the responses of a free Fermi gas in the non-Pauli-blocked region (q > 2/c_f) 
can be expressed in terms of the variable only (apart from g-dependent multiplicative 
factors). We shall see that even in the Pauli-blocked region and for an interacting system it 
is convenient to use the pair of variables (q,ip) instead of (q,uj). 

Besides the obvious advantages related to the use of a scaling variable (see Section IV), 
there is another reason for expressing the responses in terms of ip: when the latter is used the 
responses viewed as functions of ip turn out to adjust to the form assumed for the nucleon 
kinetic energy. To be more specific, starting from either a non-relativistic or relativistic 
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Fermi gas, one is always led to essentially the same dependence of the responses upon the 
corresponding ip variable. 

We shall see in the following subsections that the energy denominators of the free nucleon 
propagators appearing in the Feynman diagrams for the response functions are always given 
by uj — e^ q + \ where is the kinetic energy of a nucleon of momentum k and k < kp. 
In the non-relativistic case 



^(0)nr _ (k + qf 



(0)nr . 




(2.9) 



where 




(2.10) 



is the standard scaling variable of the non-relativistic Fermi gas and the nucleon mass. 
In the relativistic case, one would have 

" " + i° )r = « - + l) 2 + m% + ^k 2 + m%; (2.11) 

however, in Ref. [32] it was shown that at the pole (where the above vanishes) a very good 
approximation for Eq. (2.11) obtains by using Eq. (2.9) with, instead of ip nr , 



1 

kF 



u)m N {l + uj/2m N ) q 
a 2 



(2.12) 



and then by multiplying the free response by l+u/m N , which is proportional to the Jacobian 
of the transformation from the variable uj to the variable ip. Thus the use of the scaling 
variable in Eq. (2.12) entails the substitution 

(0)r (0)r ( W \ (0)nr (0)nr 



^ - + - w + - 4t; + $r- (2.i3) 

In turn, this implies that the pole (which provides the contribution to the imaginary part of 



the propagator) is located at uj — \Jm 2 N + q 2 + 2q ■ k — m N , namely at the place predicted 
by the exact expression in Eq. (2.11) when k 2 is neglected with respect to m 2 N . As stated 
above, since k is always below kp, this is a good approximation and, indeed, the free RFG 
response calculated using the scaling variable in Eq. (2.12) reproduces that of the exact 
calculation accurately, the discrepancy being typically below 1%. 

However, in the calculation of higher-order (RPA) contributions, the real part of the 
energy denominators also comes into play and the validity of the approximation far from 
the pole should also be checked. With some algebra — and assuming k 2 /m 2 N C 1 — one 
can write 
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FIG. 3. Ren(°) at q = 500 MeV/c and 1000 MeV/c as a function of the energy transfer, using 
the exact relativistic kinetic energies (solid) and the approximation discussed in the text (dashed); 
here k F = 195 MeV/c. 



(0)r (0)r ^ 

U ~ e k+q + C fc - 



qk F 



tp r — q ■ kjkp 




2 \ rriN 

qkp ip r — q • k/k F 
m N l + u/m N 



1 + 



q 2 + 2q-k 



m 



N 



(2.14) 



where, in the last passage, we have replaced the square root with its value at the pole. In 
Fig. 3, we display the real part of the free polarization propagator (defined in the following 
subsection) using the exact relativistic dispersion relation and the prescription of Eq. (2.14) 
at q — 500 MeV/c and 1 GeV/c as a function of u. The agreement between the two ways 
of calculating RelT ) is quite good at both momenta. 

Equation (2.14) provides an approximation for the free ph propagator. A prescription 
to obtain the (kinematically) relativistic polarization propagators at any order in the RPA 
expansion (see Section II E) can easily be obtained by noting that Il( n ) — the n-th order 
contribution to the RPA chain — contains n + 1 ph propagators; one then has 



n (n)r (?, 



m N 



n+l 



U^ nr (q,cu{l + cu/2m N )). 



(2.15a) 



Actually, all of the response functions derived below are expressed in terms of a generic 
scaling variable ip, as U^^q, ip). One can then get the non-relativistic response by using the 
(exact) expression in Eq. (2.10) for ijj and the relativistic response by using the (approximate) 
form in Eq. (2.12) and multiplying each polarization propagator by the appropriate power 
of 1 + cu/rriN, i. e. 
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n(">(g,cj) = (l + —Y + n (n)nr (g, i/j r ). (2.15b) 

Note that in the calculations of Refs. [14,23] only an overall Jacobian factor, 1 +uj/m N , has 
been applied to the RPA response functions. In typical kinematical conditions the size of 
the error introduced by this further approximation is of the order of a few percent. 



C. Free response 

Although the free Fermi gas response function is a subject for textbooks (see, e. g., 
Ref. [20]), it is useful to derive it here using a slightly different approach, since it illustrates at 
the simplest level the method we have adopted to overcome a technical difficulty one meets 
in nuclear matter calculations — namely the presence of ^-functions, which considerably 
complicates analytic integrations. As a side effect, the expression for Il( ) also comes out to 
be much more compact than in standard treatments. 

From Eqs. (2.5) and (2.6), one immediately finds that 

ngU = ngU = n^U = nSSU = nw (2.16) 

where following Eqs. (2.3) and (2.4) we have defined 

nW(,,a;) = /^Gg(fc,g;a;), (2.17) 

having set 

Gph {k,q;u) = -ij |^G (0) (fc + q, k + u)G<® (k,k ), (2.18) 
G^(k, k ) being the free one-body propagator 

GM(*, fa) = '<*-*'>. + ° {kF - k \ ■ (2.19) 

The integration over k in Eq. (2.18) is straightforward, yielding 

(o) A _ 9(k F -k)d(\k + g\-k F ) 9(k - k F )9(k F - \k + g[) 

U ph {K,q,u)- (Q) (o) + — — (o) (o) . > 

U ~ € k+q + e k + l7 l ~ W + e fc+q ~ 6 fc + ?7 7 

which, inserted back into Eq. (2.17), would give the standard definition of Il( ). Instead, let 
us rewrite G^ as 

(0) A _ 6(k F - k)6(\k + q\ - k F ) 6{k-k F )9{k F -\k + q\) 

%h <?> - ( o) (0) . + (0) (0) . . 

, fl(A; F - A;)#(A; F - |fc + g|) + 0(/cf — k)6(k F — |fc + gj) 2 



(0) (0) . • 1 . (0) (0) 
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having added and subtracted the quantity in the second line, where we have set r] w = 
sign(u)ri. A few algebraic manipulations then yield 



Hence, from Eq. (2.17) one gets 



9{k F - k) - 0{k F - \k + q\) 



uj — e 



(0) 
k+q 



(2.22) 



(2nf 

m N k\ 
T(2^ 



(o) (o) . • 



—uj — e 



(0) 
k+q 



(0) 



(2.23) 



Note that only one ^-function forcing k below k F is left, Pauli blocking being enforced by 
cancellations between the energy denominators. In Eq. (2.23), we have introduced q = q/k F 
and the dimensionless function 



e <0, W = i/>_4 

which is easily evaluated, yielding 



(2.24) 



ReQ (0) (^) 



^ + -(l-^ 2 )ln 



ImQ^(^) = -sign(w)0(l - ^)'-(\ - ^) = -fflgn(w)0(l - ^ 2 )^[P W - P2 



7T 



7T , 



(2.25a) 
(2.25b) 



where P n and Q n are Legendre polynomials and Legendre functions of the second kind, 
respectively. 

The expression in Eq. (2.23) has a simple physical interpretation. If one switches off Pauli 
blocking, the response of a Fermi sphere, with four particles per momentum state up to k F , 
is given by a parabola over the response region q 2 jlm^ — qk F jm^ < uo < q 2 /2mN + qk F /mN, 
that is, the curve obtained joining the dotted line and the parabolic section of the solid line 
in Fig. 4. With respect to the Pauli blocking, two kinds of spurious terms arise when k and 
\k + q\ are both below the Fermi surface. If \k + q\ > k, then a spurious contribution occurs 
in the Pauli- forbidden region < uj < qk F /m^ — q 2 /2m^, whereas if \k + q\ < k, then a 
contribution occurs with the same strength for q 2 /2m^ — qkp/m^ < uj < 0. Hence, in order 
to get the correct response function, one simply subtracts — for a given uj > in the Pauli- 
forbidden region — the total of the spurious contributions at —u, thus getting the familiar 
linear dependence on uj. Graphically, as illustrated in Fig. 4, this amounts to reflecting the 
response at negative transferred energies in the vertical axis and then subtracting it. 



D. Hartree-Fock response 

The HF polarization propagator in nuclear matter is obtained by dressing the one-body 
propagators appearing in 11^ with the first-order self-energy S^, so that one can follow 
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w (MeV) 

FIG. 4. Free response function at q = 300 MeV/c and Jcf = 195 MeV/c. The parabola given 
by the dotted line plus the parabolic part of the solid line represents the response of a Fermi sphere 
without Pauli blocking; the dashed line represents the response at negative energies reflected in 
the vertical axis and, once subtracted from the dotted line, yields the solid straight line, that is, 
the Pauli-blocked part of the response. 



essentially the same derivation of the previous subsection. The spin-isospin matrix elements 
are the same as for the free response, yielding 

LL L;I=Q — ii L;/=l — ii T;/=0 — iJ T;J=l — 11 > {Z.ZO) 

where 

nHF(? ' w)= /w^ (fc,<,;w) (2 - 27) 



and 

<-°° dk 



^ ^ GHF ( fc + 9,ko + Lu)G RF (k, k ), (2.28) 
G RF (k, k ) being the HF one-body propagator 



G»(t, U) = + , (2.29) 

with 



e^e^+E^*)- (2.30) 
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Since the first-order self-energy does not depend on the energy, the integration over ko 
can be performed along the lines of Eqs. (2.20)-(2.22), yielding 



and, finally, 

n HF (^) = /^(^-*) 



6(k F - k) -9(k F - \k + q\) 



uj 



uj — e 



(i) 

k+q 



t k+q + t k + lr lfJ 



+ 



(2.31) 



,(1) 



(1) 



+ €),' + vq u -uj - e K k ' +€),'- ir], 



[2.32) 



The HF response function is proportional to the imaginary part of II 



HF. 



m N 



ImH HF (g, u) = -sign(c)7r J j^0{k F - k) [5 (uj - e^ q + e 
= — sign(w)7r 

X / 



5 f — uj — e 



dk k 



q (2tt) 2 
1 



(i) , W\\ 

k+q + e fc J 

(2.33) 



m N 



^ k 2 + q 2 + 2qyo ) - m* N (^/k 2 + q 2 + 2qy Q ) 



having defined the effective mass as 



m N 



1 + 



m N 



(2.34a) 



k dk 



or 



rn 



N 



+ k 2 



1 + 



m 



% + k 2 dE« 



(2.34b) 



k dk 

for the non-relativistic or relativistic case, respectively, whereas y$ and y solve the equations 



/hfH^Z/o) = 
fuF{-u\k,y ) = 0, 



with 



f£?(u\k,y) = uj 



qy 



S (D(^2 + ?2 + 2?y)+E (l) W 

fn F (u\k, y) = uj- \/m 2 N + k 2 + q 2 + 2qy + ^Jm 2 N + k 2 



qy 



m N 



}-WUk 2 + q 2 + 2qy) + E^(k). 



[2.35) 



(2.36a) 



(2.36b) 



Although the evaluation of the HF response is numerically quite straightforward, in Ref. [14] 
an analytic approximation for ImIl HF has been worked out, with the aim of using it to include 
the HF field in RPA calculations. Here, it will be shown that the analytic approximation 
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is valid not only for the HF response, but more generally, although in the HF case one can 
directly assess the good accuracy of the procedure. 

In any Feynman diagram considered here and in the following, the nucleon self-energy 
enters through the ph energy denominators, 



LO 



,(l)nr (l)nr 



"fc+q 



= oo 



+ 



k 2 



(k + q) 



2,771 j\r 2777 jv 



^(\k + q\) + ^(k), 



(2.37) 



where the non-relativistic expression for the nucleon kinetic energy has been used. In 
Eq. (2.37), one can always assume that k < k F and \k + q\ > k F . Although the latter 
inequality is not immediately apparent from, e. g., Eq. (2.32), remember that cancellations 
between the energy denominators are such as to enforce the Pauli principle; the same will 
also be true for the RPA diagrams 3 . 

Clearly, if E^(fc) were parabolic in the momentum, the inclusion of the self-energy 
would be achieved simply by substituting an effective mass for m^. For realistic potentials, 
a parabolic fit for the self-energy over the whole range of momenta is in general not a good 
approximation. It is a good approximation, on the other hand, to fit separately the particle 
and hole parts of the self-energy, the fit being restricted to the range of momenta actually 
involved in the integration. Since in Eq. (2.32) (but also in the RPA diagrams discussed 
later) k is integrated from to k F and, furthermore, \k + q\ > k F , one can set 



k 2 

E^^A + B , 0<A;<A; F , 

2m N 



k 2 

E (1) =A + B , max(g - k F ,k F ) < k < q + k F . 

2m?q 



(2.38) 



Inserting this "biparabolic approximation" back into Eq. (2.37), and setting e = A — A and 
m*™ = m N /(l + B), one gets 



10 ~ e h+q + C fc 



u-(l + B)-^ (l + B) 



q ■ k 



k 2 



2m A? 



qk F 

m*^ ) k, 



777™ r q 



+ A-A + (B-B)- 

777JV 2777JV 

- A B-B k F f kV 
q '~kF~ + 1 + B2q~ [k^J 



qk F 

777^ r 



k 

^nr-q-jr 

k F 



(2.39) 



To go from the second to the last line in Eq. (2.39), we have neglected the term proportional 
to k 2 , which is expected to be small, since k < k F and, typically, q > k F . However, this 
approximation depends upon the interaction and one should check its validity, since it affects 
both the parameters B and B. In Ref. [14] the term neglected has been shown to be small for 



3 It should also be noted that the infinite Fermi gas is more in touch with the physics for relatively 
large momenta (q > 2k F ), where the above conditions are satisfied by definition. 
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the Bonn potential; the same turns out to be true also for the effective interaction employed 
in the next section. 

Equation (2.39) is similar to the expression (2.9) for the free energy denominator, but 
for the substitutions 



*nr m N 

m N -> m N = 



1 + B 



(u + e) 



rrr N q 



(2.40) 



1 + B ' 



X 



k F \ q 2 



(or uj — > + e). 

In Ref. [14] relativistic kinematics had been accounted for by applying to the above 
formulae the substitution uj — > uj(1 + ui/2m N ) previously discussed. The correct ap- 
proximation can be worked out by starting again from the ph propagator by defining 
(AE^fc, q) = E^(k) - E^flfc + q\)) and rewriting it as 



U ~ € k+q + € k 



UJ 



+ ^/k 2 + m 2 N + A^\k,q) + yf{k + q) 2 + m 



cu 2 + 2u^k 2 + mfv + 2(cu + ^A; 2 + m^AEW (k, q) + [AEt 1 ) (Jfe, q)] 2 - q 2 - 2q ■ k 
m*N 1 + u/m N + A^/m N 



qki 



ip* - q ■ k/k F 



(2.41) 



where 



m 



m N 



' 1 + B(l + u/m N + AW/mjv) 

V»p + X[l + B(l + cu/m^v + A«/2mAr)] 



A« = £ 



1 + 5(1 

,.2 



fcu/mjv + A( 1 )/m JV ) 



(2.42) 



2m 



A? 



with x already defined in Eq. (2.40). In deriving Eq. (2.41), we have assumed that k 2 <C m%, 
have evaluated the numerator at the pole thus discarding any angular dependence and, in 
the denominator, have retained only terms at most linear in q-k/kp. As one can see, besides 
the transformation uj — > u(l + u/2m N ) there are other relativistic corrections, both to the 
effective scaling variable and to the Jacobian. 

The quality of the approximations introduced above is good: indeed the HF response is 
reproduced with at most a few percent discrepancy (except on the borders of the response 
region, where the Fermi gas is anyway unrealistic). Thus, we see that in either the non- 
relativistic or relativistic case, the prescription to include HF correlations in a response 
function is simply to replace ip with ip* and m N with m* N (and to multiply by a normalization 
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factor when employing relativistic kinematics (see Eqs. (2.15)). For instance, from Eq. (2.23) 
one gets 



n HF (g,^) = J 



[2A3) 



with 



and 



Jnr 1 



r , " A(1) 
J r = 1 + + 

rriN rriN 



(2.44) 



(2.45) 



In Appendix B, we give the explicit expressions for the first-order self-energy, based on 
the generic potential in Eqs. (2.1)-(2.2). 



E. Random phase approximation response 

If in Eq. (2.3) one substitutes the irreducible vertex function T 13 with the matrix elements 
of the bare potential, one gets the so-called random phase approximation to G ph . In terms 
of the polarization propagator in Eq. (2.4) one would get an infinite sum of diagrams such 
as those shown in Fig. 5. 

We have already noted at the beginning of Section II A that, while for the two-body 
Green's function G ph one can introduce an integral equation, this is not in general possible 
for the polarization propagator. It becomes possible when one approximates the irreducible 
vertex function T 13 with the direct matrix elements of the interaction. In that case, in an 
infinite system one gets a simple algebraic equation whose solution, for the polarization 
propagators in Eq. (2.5) and the interaction in Eq. (2.1), is readily found to be 

Hx (g ' U) = 1 ttW^ ? (2 ' 46) 

1 - n x (?iWj/ni )(g,w) 
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where represents the first-order direct polarization propagator: 

ng d =0(1) (g,c) = rt \q,Lu)W 0(T) (q)rt \q,tu), (2.47a) 

n T;/ d =o(i)(9^) = rt°\q,u;nV a{aT) (q) - V t(tr) (q)]rf°\q,u). (2.47b) 

The effect of the exchange diagrams is often included through an effective zero-range inter- 
action, calculated by taking the limit q — > of the first-order exchange contribution and 
rewriting it as an effective first-order direct term [33]. Exact calculations, however, show 
that extrapolating this approximation to finite transferred momenta is not always reliable 
[6]. 

A more advanced approximation scheme is given by the continued fraction (CF) expan- 
sion [34,4,5,35]. At infinite order the CF expansion exactly corresponds to the summation 
of the perturbative series, so that it is not any easier to calculate than the exact expression. 
However, when truncated at finite order, not only does it reproduce the standard pertur- 
bative series at the same order, but in addition it yields an estimate for each one of the 
infinite number of higher-order contributions. Regrettably no general methods are available 
to predict the convergence of the CF expansion, the only reliable test being to compare the 
results at successive orders. 

On the other hand, one should note that for zero-range forces the first-order CF expansion 
already gives the exact (albeit trivial) result, making one hope that the short-range nature 
of the nuclear interactions allows for a fast convergence. Indeed, all available calculations 
have been performed truncating the CF expansion at first order [4,5,27,28,14,23]. Here, as 
anticipated, we shall test the convergence up to second order. 

The CF formalism for the polarization propagator is developed in Ref. [4] for the case of 
Tamm-Dancoff correlations and extended in Ref. [5] to the full RPA. Instead of following the 
rather involved formal derivation given there, here we shall briefly sketch a sort of heuristic 
derivation of the CF expansion. 

Let us assume that we want to build a CF-like expansion for the polarization propagator, 
according to the pattern 

TT(°) 

n RPA = . (2.48) 

I — A 

D 

1-C- 



1 - ... 

We have said that the CF approach at n-th order exactly corresponds to the perturbative 
series at the same order and then it approximates the higher orders. Thus, if we want to 
approximate the exact RPA propagator at first order in CF (for sake of illustration we drop 
spin-isospin indices), 



oo 



n=0 

we can rather naturally set 



-qO) ^ -q(O) 



n(°) 



(2.50) 
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In Eq. (2.50) = lT°)4VTl( ) + U^ ex is the sum of the direct and exchange first-order 
terms of RPA — since this yields the correct expression for the direct terms. With the 
approximation in Eq. (2.50) the summation is trivial, yielding 

ttRPA = 11 = (o *i 1 

CF1 l-nw/nc) i-4yn(°) -n^/nc)' K ' 

We could then add in the denominator of the above expression the exact second-order 
term, , after subtracting its approximate estimate given by the first-order CF expansion, 
[nWp/nW. We would thus obtain 

jjrpa _ n(°) 



( ' K2 i - nw/n(°) - {n( 2 )/n( ) - [nw/nw] 2 } 



~ i - 4yn(°) - n( 1 ) cx /n(°) - {n( 2 ) ex /n(°) - [nw^/nw] 2 }' 



(2.52) 



It is easily deduced from Eq. (2.52) that the third-order term is approximated as IT( 3 ) = 
n( 1 ){2n( 2 )/n(°) - [nW/n^] 2 }. Then, going ahead in a CF-style expansion we would guess 
for the exact RPA propagator the following expression: 

ttrpa = 

" 1 _ n(1) /TT(0) rr^ynw - [n^/nop " 

' n( 3 ) cx /n(°) + [n^^/n^] 3 - 2[n( 1 ) cx /n(°)][n( 2 ) cx /n(°)] 

1 n (2)ex^ n (0) _ [ n (l)ex/ n (0)]2 

(2.53) 

This is the expression that one would get from the formalism of Refs. [4,5] when the expansion 
up to third order is worked out. Note that we did not assume any specific scheme (either 
Tamm-Dancoff or RPA) in this heuristic derivation. 
Thus, following Eq. (2.46), we can write 

n(°) 

lf PA = (2.54) 



i - n« d /n(°) - r4 1)e 7n(°) - — x _ — 

where n^ d has been defined in Eq. (2.47). Clearly, a truncation at ra-th order would require 
the calculation of the exchange contributions up to that order. Exploiting Eq. (2.5) these 
can be cast in the form 

<} ex (g, u) = tr[6[n^( q , u)Ol\ = £ C^n^ n ( q , w ), (2.55a) 

nS ex (g,c) = EA 3i trfenW cx ( 9lW )oy = E^j^iXM, (2.55b) 

ij oa 

where the indices ctj run over all the spin-isospin channels and the spin-isospin factors 
are absorbed into the coefficients C^- an = cjj? l) c£ 2) ...c£ n) (see Table I). Moreover the 
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X 




C T x 


c x 


/~~t<JT 
I 


&x 




7;7 = 


1 


3 


3 


9 








L;I=l 


1 


-1 


3 


-3 








T;7 = 


1 


3 


-1 


-3 


-1 


-3 


T;I = 1 


1 


-1 


-1 


1 


-1 


1 



TABLE I. The spin-isospin coefficients C x (see text) in the longitudinal and transverse isoscalar 
and isovector channels, for the interaction in Eq. (2.1). 

"elementary" exchange contribution n^ e * n containing n interaction lines V ai ...V an , namely 4 

4tU?, <«>) = /!§••• d ^G^\K 1 )G^\K 1 + Q)V ai ( kl -k 2 )-.. 
■ ■ -V an (k n - k n+1 )G^(K n+1 )G^(K n+1 + Q) 

= / (S? • • • ^ G >^-) v ^ -*>)■•■ 

• • -K n (fc„ - fc n+1 )Gg(fc n+1 , 9 ;u;) (2.56) 

have been introduced. With the definition of given in Eq. (2.22) and by a suitable 
change of integration variables one can eliminate all of the ^-functions that contain angular 
integration variables, leaving a multiple integral with the following general structure: 

4tUg,c) = {-IT J ^e(k F - kl ) ■ ■ ■ d ^o(k F - k n+1 ) 

x [ — — ^(fei - k 2 ) ■ ■ ■ V an (k n - k n+1 ) — — 

+ J2(u ^ -u)]. (2.57) 

In Eq. (2.57), XX W — * ~ u0 ) stands for the sum of all the terms generated according to the 
following rules: 

i) Take all of the terms obtained by substituting u — > — u in one energy denominator in 
the second line of Eq. (2.57); then add the contribution obtained by performing the 
same substitution in two energy denominators and so on up to when the replacement 
uj — > — uj has been performed in all the n + 1 denominators; 

ii) Every time (uj - e kt+q + e ki + irj^y 1 is replaced with (-uj - e ki+q + e ki - irj^y 1 then 
replace fc, with — ki — q in the potential. 



4 The following formulae are valid for non-tensor interactions; the treatment of the tensor terms 
is slightly more complex and it is given in Appendix C. 
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The number of integrations can be reduced by noticing that the azimuthal angles are 
contained only in the potential functions V ai . For typical potentials this integration can 
be done analytically, hence it is convenient to introduce a new function representing the 
azimuthal integral of the potential. To this end, define the new variables: 

\k - k'\ = \jk 2 + k' 2 - 2kk'[cos9 cos 9' + sin6» sin 6' cos(y? - (p')} 



k 2 + k' - 2[yy> + Jk 2 - y 2 yk' 2 - y' 2 cos(f - f')] (2.58) 



= \J x + x' — 2\fx\fx' cos(y? — f') + (y — y') 2 , 
where y = k cos 9 and x = k 2 — y 2 . Then, one can introduce 

W a (x, y- x', y') = ^V a (k - k') = W a (x', y'; x, y) (2.59) 

JO Z7T 

and rewrite Eq. (2.57) as 



i i 

x- -—W ai (x 1 ,y 1 ;x 2 ,y 2 ) • • • W an (x n ,y n ;x n+1 ,y n+1 )- — — 

tp-yi + irju ip- y n+ i + vq u 

+ £>;->-£«;). (2.60) 



For n = 1 one has 



x — — W a {x, y; x', y')- — — 

ip-y + irju, ip-y + ir]u> 



'm N \ 2 kp 



q J (2tt) 4 



Qi 1} (o, 1>) - QL 1} (?, 1>) + (o, <A + ?) - Qi x) (-?, <A + ?y , 

(2.61) 



where 



Qi%, V) = 2 /* dy- ^— C dy'W a "(y,y';q) (2.62) 

J-i ip-y + irju, J -i y-y' + q 



and 



WZ(y, y';q) = - [ " dx-f dx' W a (x, y + q; x', y'). (2.63) 
2 Jo Z Jo 

Note that in getting to Eq. (2.62) use has been made of the Poincare-Bertrand theorem 
For the potential in Eq. (2.2) W" can be calculated analytically (see Appendix D), so 
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that the calculation of the first-order exchange contribution to the polarization propagator 
is reduced to the numerical evaluation of two-dimensional integrals for the real part and of 
one-dimensional integrals for the imaginary part. 
For n = 2 one has 



— J_ idyi -J Q d Xl Jjy 2 -J Q dx 2 Jjy,-J Q 
1 1 1 



q J (2tt 



dx 



3 



x -—W a (x 1 ,y l ;x 2 ,y 2 )- ——W a (x 2 , y 2 ; x 3 , y 3 )- 

V>-yi + irh ip-y2 + ir]u> W-V3 + irh 

( m N 



3 k 6 
| F 

q ) (2vr)6 



-2^(0, 0; V + g) + QS(0, -g; V + ?) + Q^(-?, 0;^ + q) - Q^(-q, -q\ i> + q)] , 

(2.64) 



where 



(2) _ f 1 1 f 1 '^ 1 

Qj*(<h,<b\$)= I dy-\ dxg a (x,y + q 1 ;ip + q 1 )———^Q a/ (x,y + q 2 ;'ip + q 2 ) 
j— i z Jo w y \ ^Wcd 



tp-y + ir] u 

(2.65) 



and 



g a (x,y;4>) = f 1 dy'- 1 ——W^x,y;y'), (2.66a) 

W^x,y;y') = \ f'" dx' W a (x, y; x', y'). (2.66b) 
2 jo 

For the potential in Eq. (2.2) W' a can be calculated analytically (see Appendix D) and one is 
left with the numerical integration of Eqs. (2.65) and (2.66a), so that the calculation of the 
second-order exchange contribution to the polarization propagator is effectively reduced to 
the numerical evaluation of at most three-dimensional integrals. Higher orders add a numer- 
ical two-dimensional integration for each additional interaction line, since, for a potential of 
the form in Eq. (2.2), only the azimuthal integration can be performed analytically for the 
interaction lines that do not close on the external vertices. 

Finally we recall that the nucleon propagators can be dressed by the HF field, as explained 
in Section II D, by replacing ip — > ip* and m N — > m* N , where ■?/>* and m* N have been defined 
in Eqs. (2.40) and (2.42), multiplying by the appropriate power of the normalization factor 
l + uj/m,]sr + A^/tojv when relativistic kinematics are employed (see Eqs. (2.15) and (2.45)). 



F. Effective particle-hole interaction 

In order to assess the contributions to the nuclear responses arising from the various 
approximation schemes introduced so far and to show typical results, first of all we have 
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to choose an effective interaction. This choice can be rather delicate as it may introduce 
uncontrolled uncertainties in the calculation. Here, however, we are not interested so much 
in comparisons with data, but rather with setting up working many-body schemes. For this 
purpose, we shall use the G-matrix based on the Bonn potential of Ref. [25], adapted to the 
quasielastic regime as in Ref. [37]. Although the attraction provided in the scalar-isoscalar 
channel by this interaction is definitely too strong [37], it will serve our illustrative needs. 

Two approaches to determine the effective ph interaction in the nuclear medium appear to 
be possible: one can either directly fix an effective potential by fitting some phenomenological 
properties or start with a bare nucleon-nucleon interaction and calculate the related G- 
matrix. Parameterizations of the ph interaction based upon the first procedure are generally 
only available at very low momentum transfers (in terms of Migdal-Landau parameters), and 
since we are probing relatively high momenta, we have resorted to using a G-matrix. We 
have chosen the one of Ref. [25], that, in our view, has the following appealing features: it 
is based upon a realistic boson-exchange potential; it accounts for the density dependence; 
and it includes (nonlocal) exchange contributions in the effective interaction, which are 
conveniently parameterized in terms of Yukawa functions. 

A feature related to the effective inclusion of antisymmetrization effects is particularly 
interesting in order to test a specific widely employed approximation scheme, the so-called 
ring approximation, in which the exchange diagrams of the RPA series are dropped and their 
effect mimicked by adding to the direct interaction matrix elements an effective exchange 
contribution (see, e. g., Ref. [33]). Indeed, below we shall compare calculations employing 
the fully antisymmetrized formalism developed in the previous subsections using the direct 
part of the G-matrix, to those employing the ring approximation using the antisymmetrized 
effective interaction. 

To facilitate the comparison with the original parameterization of Ref. [25] , the potential 
is given here using the standard representation of Eq. (2.1) in spin and isospin (no spin-orbit 
contribution will be considered in the following), but employing different symbols for the 
momentum space potentials (and adding the tensor contributions in the exchange channel): 

V(kf, hi, k F ) = F + F'ti ■ t 2 + Gcr 1 ■ <x 2 + G'cr 1 ■ <x 2 r x ■ r 2 

+TS 12 (q) + T'S 12 (q)r 1 • r 2 + HS 12 (Q) + H'S 12 {Q)t 1 ■ r 2 , (2.67) 

where q = ki — kf, Q = ki + kf (ki, kf being the relative momenta in the initial and final 
states, respectively) and the coefficients are density and momentum dependent. 

Before utilizing the interaction of Ref. [25] in a calculation of quasielastic responses, a 
few issues have to be addressed [37]. 

a) The density dependence of the G-matrix is given in terms of density-dependent cou- 
pling constants, which is not very useful for applications to finite nuclei. Furthermore, the 
parameterization is fitted for 0.95 fm _1 < k F < 1.36 fm _1 , which spans a range of densities 
down to roughly 1/3 of the central density. Extrapolation of this parameterization to lower 
densities (which is crucial for application to hadron scattering) gives unreasonable results. 
Thus, we have chosen to employ a linear p-dependence (V = V cx + V p p) , which is considered 
a reasonable choice (see, e. g., Ref. [38]). In Fig. 6 one can see a comparison of the two 
parameterizations for the fc^-dependence of the effective interaction. It should be noted that 
most of the contribution to the quasielastic responses comes from densities where the two 
descriptions differ by only a few percent. 
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FIG. 6. Effective interaction in the non-tensor channels as a function of fop at q = 0; linear 
(solid) and from Ref. [25] (dashed) density dependence. 
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FIG. 7. (a) Direct and exchange ph matrix elements; (b) direct pp matrix element. 
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FIG. 9. As in Fig. 8, but for the tensor channels; T and T' do not depend on the density. 
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b) In order to obtain a local interaction at a fixed density, one can use the relation between 
q, Q and ki, i. e., Q = \J Ak 2 — q 2 and take for ki a suitably chosen average value, (ki); then, 
the only independent momentum is q. The authors of Ref. [25] were interested in a potential 
for low excitation energy nuclear structure calculations and hence they assumed that the 
two nucleons in the initial state lie on the Fermi surface and so averaged over the relative 
angle, getting (ki) 0.7/cf- Clearly, in this case one has the constraint < q < 1.4fcp- 
On the other hand, we are interested in the ph interaction in the quasielastic region where 
one nucleon in the initial state is below the Fermi sea, while the other can be well above 
it. A look at Fig. 7 shows that ki is defined in terms of the particle and hole momenta as 
ki = (p — h')/2 = (h — h! + q)/2. Thus, at fixed q one should average ki over h and h', 
getting (ki) & \jQk 2 F /h + q 2 /2. Now ki grows with q, so that there are no longer constraints 

on q and the exchange momentum turns out to be constant, Q = ^6/5kp. In Fig. 8 one 
can see the resulting interaction in the non-tensor channels. 

c) The tensor channels are simpler, since in the parameterization of Ref. [25] there is 
no explicit density dependence (Fig. 9). The coefficients of the exchange tensor operator, 
H and H', display a very mild density dependence induced by Q, which is completely 
negligible. The only drawback concerns the treatment of Su(Q): assuming that q and Q 
are orthogonal, with some algebra one can show that Si 2 (Q) = —Si 2 (q)/2. 



G. Testing the model 

First of all we have to choose the Fermi momentum. Of course, one could easily perform 
a local density calculation to achieve a better description of finite nuclei. Here, for sake of 
illustration, we prefer to use the pure Fermi gas. The choice of kp can be made in several 
ways — here we shall choose an average value according to the formula (see, e. g., [39]) 

k F = jJ drk F (r)p(r), (2.68) 

where p(r) is the empirical Fermi density distribution normalized to the number of nucleons 
and fcfr(r) = [(3n/2)p(r)] 1 ^ 3 . For 12 C one gets kp ~ 195 MeV/c and this is the value used 
in the calculations that follow. 

Let us start with the HF response. In Fig. 10 we display the HF response of 12 C at q — 
300, 500 and 1000 MeV/c. As anticipated, Eq. (2.43) turns out to be a good approximation 
to the exact expression (2.33) (except on the borders of the response region, where the Fermi 
gas is anyway unrealistic). The HF correlations widen the response region and quench and 
harden the position of the quasielastic peak, as is well known. Note however that the short- 
range correlations, which are embodied in the effective interaction based on a G-matrix, 
reduce the amount of hardening that is observed in calculations based on the bare Bonn 
potential [14]. Note also that the same level of accuracy is obtained using either non- 
relativistic or relativistic kinematics. 

Before discussing the RPA results, we would like to test the convergence of the CF 
expansion. For this purpose, in Fig. 11 we compare the longitudinal RPA responses at 
first and second order in the CF expansion using a model one-boson-exchange interaction, 
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FIG. 10. Fermi gas responses for k F = 195 MeV/c at q = 300, 500 and 1000 MeV/c: Free 
(dotted), exact HF (solid) and HF approximated according to the prescription of Section II D 
(dashed). The ph interaction is the G-matrix discussed in the text; the kinematics are relativistic. 

V a (k) = ct\ • (T2g[m 2 / (m 2 + A; 2 )] (the spin operators having the purpose of killing the direct 
(ring) contribution). For values of the coupling constant g and of the boson mass m typical of 
realistic nucleon-nucleon potentials one finds that the first- and second-order results match at 
the level of a few percent (in the left and middle panels of Fig. 11, the solid and dashed curves 
are actually indistinguishable). One has to go to very low boson masses (a few MeV) and, 
consequently to very high values of g in order to find some discrepancies. To understand these 
results better, in Fig. 12 we display the modulus of the polarization propagator at first order, 
(dotted), at second order, Il( 2 ) (dashed) and the approximation to n( 2 ) generated by the 
first-order CF expansion (see Section HE), n( 2 ) appr = /H^ (solid). From inspection 
of the curves, it is clear that the first important element to guarantee good convergence 
is the range of the interaction. Indeed, for m = 800 MeV (short-range), H~( 2 ) and n( 2 ) a PP r 
practically coincide independent of the strength of the interaction. This, of course, should 
be expected, since for zero-range interactions the first-order CF expansion gives the exact 
result. For masses of the order of the pion mass one starts finding discrepancies between 
n( 2 ) and n( 2 ) a PP r . However, for realistic values of the interaction strength the second-order 
contribution turns out to be one order of magnitude smaller than the first-order one, and 
thus these discrepancies have little effect on the full response functions (Fig. 11). 

To understand these results it may be useful to compare the strength of the interactions 
employed here to that of one-pion-exchange, g^m 2 = / 2 /3 = 0.33 (in natural units). With 
the same units, the cases with m = 100 MeV correspond to gm 2 = 0.26 and 0.65; those 
with m = 800 MeV to gm 2 = 16.7 and 41.7; for m = 1 and 10 MeV one has gm 2 = 1.3 and 
0.65, respectively. 

To summarize, from the left and middle panels of Fig. 12 one can understand that the 
validity of the CF expansion originates from the interplay between range and strength of the 
interaction. For short-range potentials where the conventional perturbative expansion may 
not converge, the CF technique yields a good approximation for the propagators at all orders; 
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FIG. 11. Fermi gas longitudinal responses for &p = 195 MeV/c at g = 300 MeV/c, with a 
spin-spin one-boson-exchange interaction, for various values of the coupling constant and of the 
boson mass: Free response (dotted), RPA with the first-order CF expansion (dashed) and RPA 
with the second-order CF expansion (solid). Note that in the left and middle panels the dashed 
and solid lines are indistinguishable. The kinematics are non-relativistic. 
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FIG. 12. Modulus of the longitudinal polarization propagator for &p = 195 MeV/c at q = 300 
MeV/c, with a spin-spin one-boson-exchange interaction, for various values of the coupling con- 
stant and of the boson mass: First-order, IlW (dotted); exact second-order, II^ 2 ) (dashed); CF 
approximation to second order, n( 2 ) appr = IlW /U^ (solid). 
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for long-range (on the nuclear scale) forces, the CF approximation is less accurate, but the 
relative weakness of the interaction already guarantees the convergence of the conventional 
perturbative expansion. One has to go to unreasonably low masses to find a situation where 
the interaction range is very long and U.W and n^ 2 ) are of the same order (right panels in 
Fig. 12). 

We can thus conclude that the calculations of nuclear response functions in the anti- 
symmetrized RPA performed at first order in the CF expansion are indeed quite accurate. 
The same conclusion is also supported by calculations with a realistic effective interaction 
- such as the G-matrix parameterization discussed above — and including HF and rela- 
tivists kinematical effects. In fact, in Fig. 13 we show the RPA and BHF-RPA longitudinal 
responses of 12 C at q — 300, 500 and 1000 MeV/c, using the full G-matrix introduced at 
the beginning of this section. Also for the full interaction, the discrepancies between the 
first- and second-order CF responses are too small to be displayed. They are at the level of 
fractions of percent everywhere, except for the case of the isoscalar channel at 300 MeV/c, 
where they rise to a few percent due to the closeness of a singularity in the propagator 
induced by the strongly attractive interaction. Indeed, as already mentioned, the scalar- 
isoscalar channel is (too) attractive 5 and softens the quasielastic peak; the scalar-isovector 
one is repulsive and gives rise to a hardening. The effect of the HF correlations is the same 
as in the discussion of Fig. 10. In Fig. 14 the transverse response is displayed for the same 
conditions. 

Finally, it is interesting and important to test the validity of the ring approximation - 
where exchange diagrams are not included — since this approximation has been widely used 
in the literature because of its simplicity. In this scheme, the effect of antisymmetrization is 
simulated by adding to the direct interaction matrix elements an effective exchange contri- 
bution (see, e. g., Ref. [33]). For details see also Ref. [37], where a prescription to determine 
the effective exchange momentum designed for use in the quasifree region has been given. 

In Fig. 15 we display the ring and RPA responses of 12 C at q — 500 MeV/c, using the G- 
matrix parameterization. It is apparent that the only channel where the ring approximation 
works reasonably well is the spin-isovector one, which, incidentally, is the dominant one in 
(e,e') magnetic scattering; it is less accurate in all other channels, especially in the scalar- 
isoscalar one. The same considerations also apply when the HF mean field is included in the 
ring and RPA responses. Note that these results confirm those of Ref. [6], where a comparison 
of ring and RPA calculations had been done using a numerically rather involved finite nucleus 
formalism. Also in that calculation the G-matrix of Ref. [25] had been employed. 



5 As indicated by the energy position of the breathing modes; in other words, nuclear matter with 
such an interaction becomes unstable. 
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FIG. 13. Fermi gas longitudinal responses for k,F = 195 MeV/c at q = 300, 500 and 1000 
MeV/c, with the G-matrix discussed in the text: Free response (dotted), RPA (dashed) and 
BHF-RPA (solid). The kinematics are relativistic. 
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FIG. 14. The same as in Fig. 13, but for the transverse channel. 
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FIG. 15. Fermi gas responses for &p = 195 MeV/c at q = 500 MeV/c, with the G-matrix 
parameterization discussed in the text: Free response (dotted), ring approximation (dashed) and 
RPA (solid). The kinematics are relativistic. 



III. PARITY- VIOLATING ELECTRON SCATTERING AND AXIAL 

RESPONSES 

A. The asymmetry, the currents and the RFG responses 

A new window on the inclusive nuclear responses that allows us to unravel aspects of 
nuclear and nucleon structure that are otherwise inaccessible to unpolarized probes is offered 
by parity-violating electron scattering from nuclei. See Ref. [40] for a general review of the 
subject. Experiments of this type exploit longitudinally polarized electrons to measure the 
helicity asymmetry A, defined as the difference between the inclusive nuclear scattering of 
right- and left-handed electrons divided by their sum, namely 

d 2 a + — d 2 a~ 
d z cr + + d l o 

A arises from the interference between the electromagnetic current that is purely vector 
(V) and the weak neutral current that has both vector and axial-vector (A) components. 
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FIG. 16. Single-photon-exchange and Z°-exchange diagrams for polarized electron scattering 
from nuclei. 



Diagrams for the associated amplitudes in leading order of the bosons exchanged (the photon 
and the Z°) are displayed in Fig. 16. In this approximation Eq. (3.1) can be cast in the form 

_ . v L R£ V (g, v) + VTR$ V (q, u) + v t >I%Hq, u) \ 

•A- — TT~, \ T\ — 7 \ • \y-£) 

v L R L {q,u;) + v T R T {q,u) 
The numerator is parity- violating (PV) and the denominator parity-conserving (PC). Here 

-=(f) 2 (3 - 3) 



Q 2 



+ tan 2 - 



(3.4) 



and 



9 

vt> = tan 

2 \ 



91 

r,2 



+ tan 2 - 



(3.5) 



are the usual leptonic kinematical factors, 6 is the electron scattering angle and, as before, 
Q 2 = id 2 — q 2 is the spacelike four-momentum transferred from the electron to the nucleus. 

In Eq. (3.1) the nuclear (and nucleon's) structure are embedded in the electromagnetic 
PC nuclear responses R^ and R^ discussed in the previous sections, while their PV analogs 
R^ v , R^ v and are discussed in this section. Here the first (second) index in the 
superscript refers to the vector (V) or axial (A) nature of the leptonic (hadronic) WNC. 
For brevity we shall often simply refer to these by their hadronic character, i.e., the L and 
T PV responses are called "vector" and the T' response "axial". Finally the scale of the 
asymmetry is set by the factor 



A 



V2G F m 2 N \Q 2 



Am 2 



6.5 x 1(T 4 t, 



(3.6) 



N 
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which is defined in terms of the EM (a) and Fermi (Gj?) coupling constants. If there were 
no additional dependence on q and u>, then the expression for Ao would imply that the 
asymmetry grows with r = \Q 2 \/Am 2 N and hence it is not surprising that the first parity 
violation in electron scattering was observed at high energies at SLAC [41,42]. On the other 
hand, it is also clear that selective processes such as elastic scattering do contain additional 
dependences on (q, u) via form factors that may make measurements at large r extremely 
difficult. In fact, only a very few have been performed to date. Even more challenging, 
but not impossible, are experiments whose goal is to disentangle in Eq. (3.2) the separate 
contributions of the PV responses R^ v , R^ v and i?^, A . 

In the investigation of the PV nuclear responses — these can assume positive as well 
as negative values — of central importance is the isospin decomposition of the hadronic 
four-current 

(4 A ^ (0) (4,o + ^(^)i 1 o (3-7) 

into isoscalar (/=0) and isovector (/=1) components (/z is the Lorentz index). In the Stan- 
dard Model at tree level the coefficients in Eq. (3.7) in the EM sector read 

$v]em = Pv^em = 2' (3-8) 

i.e. one has the usual responses i? L and i? T . In the WNC sector, instead, they read 

PvIvnc = -2 sin 2 6 W ~ -0.461 (3.9) 
Pv]wnc = 1-2 sin 2 9 W ~ 0.539 (3.10) 

for the vector coupling and 

(3 ( a = (3.11) 
= 1 (3.12) 

for the axial one. The above results obtain with the following value of the weak mixing angle 

sin 2 6 W = 0.23055 ± 0.00041. (3.13) 

The importance of isospin becomes clear if one assumes it to be an exact symmetry in 
nuclei (which is of course only approximately true). Then for PV elastic scattering on spin 
zero, isospin zero nuclei, Eq. (3.2) reduces to 

A 

— = 2sin 2 w (3.14) 
Ao 

(see Eq. (3.25) for the coefficient of the axial leptonic current), which suggests using PV 
experiments as a tool for testing the Standard Model in the low-energy regime (see, for 
example, Ref. [40]). 

To this point we have been assuming that the strangeness content in the nucleon or 
nucleus is negligible. If this is not the case, then it has been realized that PV electron 
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scattering can be invaluable in exploring this aspect of hadronic structure (again, see Ref. [40] 
for a review that contains discussion of this issue). Indeed, when strangeness content is taken 
into account Eq. (3.14) is modified as follows 

^ = 2sin 2 ^ + -fy. (3.15) 

In the above G$ = G Ep +GE n and G$ are the electric isoscalar EM and strange form factors 
of the nucleon (the indices p and n refer to the proton and the neutron, respectively). For- 
mula (3.15) will be exploited to extract G%^ in an experiment planned at CEBAF involving 
elastic PV scattering from 4 He. 

A further clue to the strangeness content may be seen in studying PV elastic scattering 
from the proton: in this case Eq. (3.2) becomes 

A eG Ep G Ep + tG Mv Gm v + 5G Mp G Ap Ci-\fi\ 
A " eGl + rGL ' [6 ' Dj 



where 



and 



1 + 2(1 + r) tan 2 (#/2) 



(3.17) 



5 = (1 - 4 sin 2 e w )y/r(l + r)(l- e 2 ) . (3.18) 

In Eq. (3.16) G Ep , Gm p and G Ap are the electric, magnetic and axial weak form factors of 
the proton, that read 



G 



E„,M V 



2 



,WNC ,wnc\ Ge p ,m p ,WNC ,WNC 



j(0) n{l) 



G 



E„,Mn 



±(l-4sin 2 6 w )G EptMp - l -G EnMn - l -G% (3.19) 



and 



G Ap = \G Ap - l -G An - l -G { t (3.20) 

Equation (3.16) has already been exploited in experiments performed at Bates (SAM- 
PLE) and CEBAF (HAPPEX) (see Refs. [43] and [44], respectively) for unraveling the 
strangeness content of the proton. Figure 17 may help in understanding the magnitudes of 
the quantities given above. There one infers that at backward scattering angles (SAMPLE) 
it is mainly the magnetic strangeness that is measured, however with some contamination 
arising from the axial contribution, whereas at forward angles (HAPPEX) it is a mixture 
of electric and magnetic strangeness that is observed, the axial contribution being totally 
negligible in that case. 

Although the above discussions of strangeness content are focused on the responses of 
the nucleon rather than of the nucleus, in fact it is also relevant for the latter, since the 
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FIG. 17. The kinematical coefficients r, e and 5 versus the scattering angle 9; the energy of the 
incoming electron is 3.355 GeV. 



nuclear responses measured in the quasielastic regime are indeed affected by the nucleon's 
isoscalar form factors and these in turn are affected by strangeness. In fact this is an excellent 
example of how the theoretical predictions in nuclear many-body and particle physics are 
interrelated. 

In particular, with regard to the responses of the nucleus, PV experiments offer the 
opportunity of 

i) disentangling the isoscalar and isovector contributions to R^ v , R^ v and i?^ 4 , 

ii) exploring the Coulomb sum rule separately in the isoscalar and isovector channels (see 
Section IV), 

iii) measuring the neutron distribution in nuclei (see Ref. [45]; not discussed here), 

iv) investigating the nuclear axial response especially in the A region (also not discussed 
in the present work) 

and 

v) unambiguously revealing the role of the pion in nuclear excitations through the (pos- 
sible) existence of a zero in the frequency behavior of the asymmetry. 

To see how this occurs let us first split the PC and PV responses into their isospin 
components (we are neglecting the strangeness content at this point) according to 
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Rh,T — Pv^EmRI,T + /^V,EAf-^L,T (3-21) 
dAV _ n(0) r i=q Ml) d/=1 (a no) 

for the vector channel. The PV axial channel is purely isovector at tree level. Next let us 
focus on the RFG model where the longitudinal isoscalar response is essentially proportional 

to (pP^ = (Ge p + G En ) and the isovector one to (G^e) — (pE p — Ge„) ■ Since G En is 
small, especially at low r, it follows that 

RF-Rt 1 . (3.23) 

But then it is clear that the RFG PV longitudinal response is almost vanishing because of 
the opposite sign and approximately equal magnitude of the coefficients in Eqs. (3.9) and 
(3.10). This dramatic consequence of the Standard Model is displayed in Fig. 18, where the 
five responses entering in the definition of the asymmetry in Eq. (3.2) are shown for |q|=300, 
500 and 2000 MeV/c together with A. Beyond the fact that R£ v is very small, one also 
observes in the figure that 

i) a similar cancellation does not occur in the transverse channel, where in fact the 
isoscalar and the isovector responses are quite different because now (G$) 2 oc (/i p + 
/i„) 2 ~ 0.77 and (G$) 2 oc - /i n ) 2 ~ 22.1; 

ii) R^ v and R^ are negative, this being related to the sign of the axial coefficient of the 
leptonic WNC 

3n(k', s'; k, s) WNC = u(k', s')(a vl , + a Al , ltl )u(k, s), (3.24) 

where u(k, s) is the Dirac spinor of the electron. Indeed, according to the Standard 
Model, 

a v = -(1 - 4 sin 2 6 W ) ~ -0.092 and a A = -l] (3.25) 

iii) as a consequence the asymmetry is negative as well, reflecting the left-handed nature 
of the weak interaction; 

iv) the asymmetry, as previously mentioned, grows with r and 6. 

How do interactions among the constituents of the RFG modify the above predictions? 



B. The role of the pion and other mesons 

In answering the last question of the previous section we extend our model from a strict 
RFG framework to one where pions are also included [27], because then we can at least 
approximately preserve the two major requirements of Lorentz covariance and gauge invari- 
ance. Indeed the RFG cross section is built (apart from overall kinematical factors) from 
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FIG. 18. Relativistic Fermi gas response functions and asymmetry for 12 C at g=300 (a), 500 
(b) and 2000 (c) MeV/c shown as functions of w. The location of uo = \\Q 2 \/m^ is indicated 
by the arrow. Upper, the EM responses i?L and Rt are displayed. Middle: the parity-violating 
responses R^ v , R^ V and R^fi are shown. Lower, the asymmetry A is given for three values of the 
scattering angle 9. 
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the contraction of the leptonic and hadronic Lorentz tensors and is therefore a relativis- 
ts invariant, although the partition into longitudinal and transverse responses depends, of 
course, upon the reference system. Moreover the single-nucleon four-current entering into 
the RFG nuclear tensor is conserved and hence the non- interacting RFG is gauge invariant. 

When the nucleon-nucleon interaction carried by the pion is switched on it is not obvious 
that the two above mentioned properties are retained. Indeed, when the correlations and 
meson exchange currents (MEC) associated with the one pion exchange potential (OPEP) 
are introduced, as shown in Ref. [27] the Lorentz covariance and gauge invariance are violated 
(however, only slightly so) due to the following approximations that are usually made: 

i) the pion propagator is assumed to be static, 

ii) and a non-relativistic expansion of the two-body currents is performed [32]. 

In order to achieve a treatment of forces and currents that is as consistent as possible we 
first limit ourselves to the study of diagrams with only one pionic line, namely, we work in 
the first perturbative order in the N-N interaction. The correlation diagrams to be evaluated 
in this scheme are the so-called self-energy and exchange contributions, that when iterated to 
infinite order generate the previously discussed HF and RPA series, respectively. Note that 
the tadpole and ring diagrams vanish due to the spin and isospin structure of the OPEP. 
Concerning MEC, three contributions occur, the pion-in-flight term, the contact term and 
the one associated with the A. A direct comparison with the exact relativistic calculation 
[46] shows that the non-relativistic expansion of Ref. [32] is indeed quite accurate up to 
momentum transfers of the order of 1 GeV/c. 

The outcome of this is that sizable pionic contributions to the EM longitudinal (spin 
scalar, a — 0) and transverse (spin vector, a — 1) nuclear responses are found. In both cases 
the correlation effects produce a hardening of the responses, that is, a shift of the strength to 
higher u. The PV longitudinal and transverse correlated responses are simply obtained from 
the EM ones through the isospin rotations implied by the structure of the WNC discussed 
in the previous subsection; the axial-vector response will be treated separately in Section 
IIIC. 

The main points emerging from this analysis are: 

a) In isospace the contribution of the self-energy diagram to the charge response is al- 
most equally split between isoscalar (I = 0) and isovector (1=1) components. The 
latter, on the other hand, is of course overwhelming in the transverse response, due 
to the dominance of the isovector magnetic moment. In contrast, in the case of the 
pionic force the 1 = part of the exchange diagram turns out to be three times as 
large as the 1 = 1 one in the charge response and this imbalance, that becomes even 
stronger in higher orders of perturbation theory, is further strengthened by the dif- 
ference between the isoscalar and isovector form factors. The isoscalar dominance of 
the pionic exchange correlations has dramatic consequences for the PV longitudinal 
response function, as may be seen in Fig. 19, where this response is displayed as a 
function of uj with and without pionic correlations. 

A physical interpretation of why R^ v is small in the independent-particle model and 
why isospin-correlations are so important in determining its ultimate size follows from 
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FIG. 19. The longitudinal weak neutral current response R^ v at g=500 MeV/c, showing the 
free RFG (dashed) and correlated (solid) results. 

the expression for the observables in a language that explicitly refers to neutrons 
and protons rather than employing isospin labeling. Indeed, by inspecting Fig. 20, 
where the diagrams describing both the EM and PV longitudinal responses for a free 
system are displayed, one easily understands why in the non-interacting case the EM 
longitudinal response turns out to be substantial: both of its vertices can in fact be 
large, as they involve the coupling of a longitudinal photon to a proton. In contrast, 
one of the vertices entering in the non-interacting PV response is always small, since 
either the longitudinal coupling of a photon to a neutron or of a Z° to a proton is 
involved. This last fact is often phrased by saying that the 7 is blind to neutrons and 
the Z° is blind to protons (i.e., in the longitudinal channel). To quantify the meaning 
of "large" and "small" we note that typically the 7 — n coupling is about 1/10 that of 
the 7 — p, and likewise the Z° — p coupling is about 1/10 that of the Z° — n. 

The exchange correlations corresponding to the exchange of an isovector charged meson 
between the particle and hole convert a neutron (proton) into a proton (neutron), and 
thus give rise to a diagram where both couplings are large — hence the crucial role of 
such isovector correlations in determining R^ V ■ 

The above arguments clearly do not apply to the transverse response both because in 
this case the 1 = channel is much weaker than the 1 = 1 channel, being essentially 
proportional to the squares of the very different isoscalar and isovector magnetic mo- 
ments of the nucleon, and furthermore because protons and neutrons can both couple 
strongly to photons via their (comparable) magnetic moments. Being also essentially 
isovector, the axial-vector response likewise does not display the sensitivity expected 
for the longitudinal response. 

b) While the tensor component of the OPEP never contributes to the self-energy in a 
translationally invariant system, the exchange diagram gets a tensor contribution, but 
only in the transverse channel and mostly via the backward-going graphs. This implies 
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FIG. 20. Feynman diagrams representing the free particle-hole polarization propagator for 
the EM (a) and PV (b) longitudinal response. The excitation of proton (p) and neutron (n) 
particle-hole pairs is shown separately. The labels strong and weak refer to the strength of the 
nucleon coupling to the photon 7 or to the vector boson Z°. 



a different role for the pionic force in the two EM responses, a finding that should be 
tested against experiment. 

c) When we extend our analysis from first to infinite order of perturbation theory, thus 
generating the HF and RPA responses, the results do not change substantially. Al- 
though the pionic interaction is strong, it nevertheless therefore appears that for 
quasielastic kinematics its effects are reasonably small at not too small q, thus ren- 
dering perturbation theory quite accurate already at the lowest order, at least for the 
classes of diagrams studied here. 

d) The contribution of the central part of the pionic interaction to the exchange diagram 
stems largely from the 5-force and not from the finite-range one. Notably for pointlike 
nucleons the continuity equation is obeyed by the OPEP and by the related pionic 
MEC, but is not affected by its 8 component. However in keeping with the usual 
approach taken in studies of pionic effects, we include a 7rNN vertex function T n , 
whose scale is set by a mass parameter A^, to smear out the 5-piece of OPEP. As a 
consequence the continuity equation is modified by the presence of and to restore 
its validity additional MEC should be introduced. Significantly these counterterms 
affect only the pion-in-flight current, which is tiny, and therefore are quantitatively 
negligible. 

e) In contrast to the exchange diagram, the self-energy term gets a contribution only 
from momentum-dependent forces, and therefore an unmodified pionic (^-interaction 
does not contribute in this channel. In fact, the contribution coming from the 8- 
function in OPEP, when modified by the 7rNN vertex form factor, contains an effective 
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momentum-dependence and so is nonzero. 

f) Most remarkably when the interactions carried by heavier mesons (namely p, a and 
uj) are switched on, the previous conclusions on the PV responses (in particular the 
dramatic enhancement of R^ V ) are not substantially changed, as illustrated in Ref. [14], 
thus showing the central role played by the pion in the nuclear dynamics for quasielastic 
kinematics. This is accomplished, however, through rather subtle aspects of the nuclear 
many-body problem, where the interference between the pion and the other mesons 
turns out to be crucial. 



C. The axial response and the asymmetry 

Let us now turn to a discussion of the asymmetry in the pion-correlated Fermi gas 
model. For its evaluation the axial-vector response function R^A is needed. In a non- 
relativistic context the latter is related to the isovector component of the EM transverse 
response, R\ VI=l , through the simple formula [28] 

B& A (q,u) = a v ^-R%r\q,w) , (3.26) 
Gm k 

which can be extended into the relativistic regime via the prescription 

R& A (q,u>) = a v ^rJ^±FSf\q, U ) ; (3.27) 

these agree to better than 2% in the momentum range 300 MeV/c < q < 1 GeV/c for 
the RFC Although not immediately apparent, Eq. (3.27) has been demonstrated to be 
preserved even in presence of pionic correlations. The effect of pionic correlations and MEC 
effects in this response is a "hardening" (a shift to higher a;) of the peak of the response at 
intermediate values of q, which then fades away and even leads to a slight "softening" of the 
response at the highest momentum transfers considered. 

We are now in the position to calculate the asymmetry A, which is displayed in Fig. 21 
as a function of u for three values of q and for forward {9 = 10°) and backward (9 = 170°) 
scattering angles. 

Upon examining Fig. 21, we notice the significant effect occurring at moderate q (say 
300-500 MeV/c), small u and forward angles. As previously discussed, it is related to the 
large negative value assumed by the correlated R^ v , which leads to a pionic asymmetry 
that is an order-of-magnitude larger than the free RFG one. However, as u increases R^ v 
rapidly decreases until it changes sign, while R^ v stays negative; accordingly, they largely 
cancel in the numerator of the ratio expressing A and this becomes substantially lowered. 

Interestingly, an energy is reached (about 60 MeV for q = 300 MeV/c) where the cor- 
related and free RFG values of A coincide. At still larger uj a further reduction of A is 
seen to occur until at about 90 MeV it nearly vanishes. This constitutes an example of a 
dynamical restoration of a symmetry (here the left-right parity symmetry) and reflects the 
complex nature of the PV longitudinal response. The near-vanishing of the asymmetry at 
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FIG. 21. The w-dependence of „4 for different kinematical conditions. The relativistic Fermi 
gas model results are labeled RFG, while the ones with pionic effects included are labeled tt. The 
two families of curves have <j^=1.26 (solid), 1.26+10% (dashed) and 1.26-10% (dash-dotted). 
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to pa 90 MeV in Fig. 21 stems from the cancellation between the positive contribution it gets 
from R^ v and the negative one it gets from R^ v ■ This trend of course fades away at larger 
6, where the role of the longitudinal PV response gradually becomes irrelevant. Finally we 
see that at larger momenta, where the impact of correlations is no longer so strongly felt, 
the nearly perfect restoration of the left-right symmetry does not show up anymore. It is, 
however, still true (even at 1 GeV/c) that an energy exists where the free and the correlated 
values of the asymmetry coincide. 

From the lower panels of Fig. 21 it clearly appears that at backward angles A is almost 
equally sensitive to the nuclear dynamical effects and to the uncertainties in the nucleonic 
form factors. This raises the question of whether or not the asymmetry itself is a suitable 
observable for extracting information on either the nuclear correlations or the single-nucleon 
form factors. As a consequence we are led to consider three different integrated observables 
that have been introduced to emphasize one of the two aspects of the problem. They are: 

1 r fUQEP fWrnax n 

AA(q,e) = — du;A{6;q,uj)- duA{6; q,u) (3.28) 

Ao> L L mi „ Juiqep 



and 



A(q,6) = — dwA(9;q,u) (3.29) 

ALU Ju>min 



" rr;^ WEM ^ <")/*t(o, r, ^ VF ) ■ [ } 



Here w m j„ and uj max are the RFG response boundaries for a fixed q 



Umin = J (k F - q) 2 + m% - Jk% + m 2 N (3.31) 



Umax = \f (k F + q) 2 + m 2 N - \]k 2 F + m 2 N , (3.32) 
the energy interval Ao; is 

A(xJ = U) max (drain (3.33) 

= y/{k F + q) 2 + m 2 N - ^{k F - q) 2 + m 2 N (3.34) 

^ 2k F q/^jq 2 + m 2 N (3.35) 

and the hadronic functions 

W EM = v L R L + v T R T (3.36) 

W PV = v L R£ v + v T R$ v + v T ,RV? (3.37) 

are divided by 
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X T (9,T,ij;r] F ) 




(3.38) 



X T (9,T,ij]7] F ) 



a A v T 2tG m G m + 



GeGe + tGmGm 



(3.39) 



in order to extract the single- nucleoli content from the many-body one [47]. The quantity 
A is given later in Eq. (4.22): since A ~ rjp <C 1, the contributions containing A may 
often safely be neglected. Removing the single-nucleon content in this way has the following 

advantages: 

i) The pionic correlations are particularly felt by AA, as is clearly apparent from Fig. 22. 
Indeed, there we first observe that at q = 300 MeV/c the results obtained with the 
free RFG model almost vanish because of the nearly perfect cancellation between the 
contributions where u < ujqep and those where oj > ujqep', however, at larger q this 
cancellation becomes less complete, owing partly to the role played by the nucleonic 
form factors and partly to the RFG model itself, whose responses (in contrast to the 
non-relativistic case) become less and less symmetric as q increases. It is also clear that 
the correlations, in particular the exchange diagram, dramatically alter the prediction 
of the free RFG, yielding a huge AA W at small 9. This result is simply interpreted by 
observing that of the nuclear responses that enter in the asymmetry the pion has its 
greatest effect on R^ v and although in the RFG model the latter accounts only for 
at most about 10% of the total asymmetry (and this only in the forward direction), 
nevertheless the impact of the pionic correlations is violent enough to induce a large 
negative value of R^ v at small u, which is in turn reflected in the large negative value 
of AA n at small 9 displayed in Fig. 22. We deduce from this that the characteristic 
behavior of AA with 9 shown in the figure represents one of the most transparent 
signatures of pion-induced isoscalar correlations in nuclei (we recall that R\f" is purely 
isovector and that in R^ v the isoscalar contribution is strongly suppressed — see 



The observable AA has been specifically devised to enhance the signal for nuclear 
correlations and to minimize the sensitivity to the single-nucleon form factors. This 
can be inferred by observing the three (almost overlapping) lines for each family of 
curves in Fig. 22, corresponding to a variation of the strength of the effective axial- 
vector coupling, g^p, of ±10% around the canonical value = 1-26. The impact on 
AA of pionic correlations is seen to be more than an order-of-magnitude larger than 
that arising from variations in the axial-vector form factor. Similar results are found 
for the magnetic and electric strangeness form factors (see Ref. [28]). 

ii) The energy-averaged asymmetry A., on the other hand, has been devised to minimize 
the sensitivity to the pionic correlations, as these tend to cancel out in the symmetrical 
integral in Eq. (3.29). The MEC contribution, on the other hand, does not average 
out in A, although in the ph sector of the nuclear excitations it turns out to have a 
rather insignificant effect. Typical results are shown in Fig. 23 for q = 500 MeV/c. 

In particular, as seen in the expanded view in the figure, at very backward scattering 
angles where one might hope to determine the effective axial- vector coupling (again 



Ref. [27]). 
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FIG. 22. The quantity AA shown as a function of 9 for 9=300 (a), 500 (b) and 1000 MeV/c 
(c). The labeling of the curves is as in Fig. 21. 
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FIG. 23. The energy-averaged asymmetry A displayed as a function of 6 for q=50Q MeV/c, 
with the same labeling of the curves as in Fig. 21. Panel (a) shows the entire angular range, while 
(b) only shows the backward-angle region in greater detail. 



variations of ±10% around 1.26 are shown in the figure) the free RFG and pionic 
correlated results for the energy-averaged asymmetry come together (for g=500 MeV/c 
at 9 m 147°, which however varies with q). Since this special condition is presumably 
model-dependent, it is unlikely that one can count on using such particular kinematics 
to effect a determination of through the variations shown in the figure. 

In contrast to the backward-angle situation, at forward scattering angles where the 
pionic correlations induce drastic modifications in R^ v , as we have seen, here the two 
families of curves differ although certainly not as much as in the case of AA. In other 
words, the observable A has some of the properties that we are looking for when we 
construct quantities that suppress the effects of correlations while bringing out the 
dependences on the single-nucleon form factors; however, this particular observable 
appears not to be entirely optimal. Since the PV longitudinal response is so strongly 
affected by the presence of pionic correlations, it is necessary to adopt an alternative 
approach to minimize these effects, leading to the introduction of the observable in 
Eq. (3.30). 
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FIG. 24. The quantity K versus 9 for ^=300 (a), 500 (b) and 1000 MeV/c (c). The curves are 
labeled as in Fig. 21. 

iii) The procedure proposed in Refs. [28,47] for the definition of the quantity 1Z has the 
goal of scaling the results (see the next section) by dividing out most of the single- 
nucleon content through the use of the dividing factors Xt and X^. In Fig. 24 we show 
1Z as a function of 6 for three values of q. Again two families of curves are displayed, 
one for the free RFG and one for the model with pionic effects included (labeled n), 
and each family has three curves (g^ = 1-26, 1.26 + 10% and 1.26 — 10%). In panel 
b we see a significant range of angles over which the pionic effects provide negligible 
modifications with respect to the RFG results and where the (merged) curves with 
the three values of the isovector/axial- vector strength can clearly be discerned. This 
behavior is similar at q — 1 GeV/c, although not quite as nicely separated. Even 
so, the difference between the RFG and pionic families at backward scattering angles 
amounts to an effective change in of only about 4%. 

It thus appears that the observable 71 is sufficiently uncontaminated by correlation 
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effects for favorable kinematics and better suited than A to disentangling the nucleonic 
form factors from the nuclear dynamics. It has been specifically designed to integrate 
out the anomalous w-dependence, leaving quantities that for the most part only retain 
sensitivities to variations in the single-nucleon form factors. 

IV. SCALING AND SUM RULES 

In this section we briefly address several issues that arise in discussing scaling and sum 
rules, showing how these two properties are interrelated and how they constrain the nuclear 
models. At sufficiently large three-momentum transfer q the so-called y-scaling region occurs 
when the energy transfer oo is lower than its value at the quasielastic peak and when non- 
quasielastic processes such as meson production do not affect the nuclear responses. The 
y-scaling approach attempts to find some function of q and cu, here denoted G, such that 
when divided into the inclusive electron scattering cross section, the result is a reduced 
response F(q,y) that scales as a function of y. Here y is an appropriately chosen scaling 
variable (see below), is a function of (q, u) and replaces cu (i.e., one uses the variables q and 
y, rather than q and uj). "Scaling" means that for sufficiently large momentum transfers the 
function F becomes universal, namely a function only of y, but not of q: 

F(q,y) q ^F(y) = F(oo,y). (4.1) 

The choice of the dividing function and scaling variable must be such as to remove the single- 
nucleon content from the nuclear responses in as model-independent a manner as possible 
while still retaining essential relativistic effects whenever feasible. 

A parallel strategy concerns the Coulomb Sum Rule (CSR) [39,48]: a dividing function 
_£/l can be devised such that the corresponding reduced longitudinal response t*l = R^/H^ 
fulfills the CSR. 

In past work medium- and high-energy data have been tested with both the usual y- 
scaling approach (for a review, see Ref. [49]), while more recently the RFG-motivated ap- 
proach has been applied and seen to scale successfully [3,50,51]. Actually these data appear 
to support not only scaling, but also superscaling, namely the existence of a function related 
to F that is the same for all nuclei [47,50,51]. Additionally, it now appears that the experi- 
mental CSR is reasonably well saturated at high momentum transfers [1] . As a consequence, 
any reliable nuclear model should simultaneously fulfill the two major requirements of 

1) scaling 

2) fulfilling the Coulomb Sum Rule. 

We now show that the RFG model simultaneously satisfies the above properties. In the 
context of the PWIA for (e, e'N) reactions the y-scaling variable is defined to be equal and 
opposite to the smallest value of the missing momentum, p m i n , attained in the y-scaling 
region: It turns out to read 

V = -Pmin = { (M° A + U>) \jw* ~ (M%-1 + ™n) - - msf 
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- q 

where 



W 2 + (M° A _^f -m 



(4.2) 



W = J(M\ + uf -q\ (4.3) 



M A and M A _ X being the masses of the initial and daughter nuclei (in their ground states), 
respectively. The energy transfer can, of course, be expressed in terms of q and y. In 
particular, it must lie in the range [39] uj t < uj < q, where 

oo t = E s + + m N f + g 2 - (M°_ x + m N ) (4.4) 

is the threshold energy, with E s = m N + M A-1 — M A being the nuclear separation energy. 
The scaling variable in Eq. (4.2) vanishes when 

uj = ojq = E s + \J m% + q 2 - m N , (4.5) 

which is roughly the position of the quasielastic peak, and hence the scaling region is char- 
acterized by having y negative. 

In order to study the scaling behavior of the inclusive (e, e') process one should be able 
to remove the effective eN cross section from under the integrals involved in going from 
coincidence to inclusive scattering. These integrals extend over the missing momentum (p m ) 
and over an energy that characterizes the degree of excitation of the daughter nucleus, £: 

£ = E A ^ - E° A ^ > 0, (4.6) 

where Ea-i is the energy of the unobserved daughter system (in general in an excited state) 
and E A _ X is that energy when this system is in its ground state, i.e., has mass M A _ V 
Naturally £ can be re-expressed in terms of the missing energy E m ; one has roughly that 
E m = £+E s . If one assumes that the proton and neutron distributions inside the nucleus are 
equal, which is a reasonable approximation for N = Z nuclei, and that the most important 
contributions to the nuclear spectral function arise from the lowest values of (p, £) that can 
be reached for given values of q and y (in the scaling region these are £ = and p = —y), 
then the function one hopes will scale as a function of y when q — > oo is 

_ l/ . d 2 aldVL e dui , A 
F (q, y) = — ( P nV ( 4 - 7 ) 

<*eN(q,y;p= -y,£ = o) 

In the PWIA this is indeed found to be the case, since in the limit q — > oo Eq. (4.7) becomes 
a function only of y, namely it scales [49]. 

Turning now to the RFG model, its spectral function [39] is 

S RFG (P, £) = ^0(k F - p)5 [£{p) - £ RFG (p)} , (4.8) 

where the excitation energy is 
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£ RFG (P) = {^F + m 2 N - y/p 1 + m%) . (4.9) 

Defining the RFG scaling variable through the intercept of the support of the RFG spectral 
function given in Eq. (4.8) and the kinematical boundaries in the missing energy-missing 
momentum plane one obtains the ^/-scaling variable of RFG 



Urfg = ™n( = m N 1 + -j- - kJ , (4.10) 

where the dimensionless variables A = u/2rriN, k = q/ZmN and r = k 2 — A 2 have been 
introduced. 

A different scaling variable was originally proposed for the RFG in Ref . [47] , namely 

^ = 4r i X ~ T ; ( 4 - n ) 

V ^(l + A)r + /^r(l + r) 



where £^ = 6^ — l = yl+77 2 ? — 1 and r\p = kp/rriN are the dimensionless Fermi kinetic 
energy and momentum, respectively. With some algebra it can be shown that the relation 
between the two scaling variables is 



£W> 2 = V 1 + (yRF G /m N y - 1. (4.12) 

The physical significance of ijj is then immediately apparent: among the nucleons responding 
to an external probe one has the smallest kinetic energy and this is given by ip 2 (in units 
of the dimensionless Fermi kinetic energy £p). Instead of working from the unseparated 
inclusive cross section towards a reduced response that, if successful, would scale as q — > oo, 
one can work directly with the separated longitudinal and transverse responses, and Rt, 
since 

1) we are most interested in model-to- model comparisons and the same procedures may 
be followed in each case (i.e., focusing on L or T responses directly), 

2) a few cases exist where L/T separations have been performed experimentally, 

3) we wish to draw comparisons with studies of the CSR where only the L response is 
relevant. 

In this spirit we seek reduced responses denoted F L)T («, ^) that scale. These are to be 
obtained from the inclusive response functions /?l,t( k , A) by dividing through by specific 
functions, denoted Gl,t(«, A): 

F LlT («, V) = iV(«> A)/G L) t(k, A )- ( 4 - 13 ) 

If the dividing functions are chosen appropriately, then as above the reduced responses 
defined in Eq. (4.13) will scale, namely, become functions only of a single scaling variable 
such as ip defined above when k — > oo, 
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Such dividing functions, derived in Ref. [52,39], are 

ZU Lp + NU Ln 



G l (k,\) 



and 



G T («,A) 



where (see Ref. [53]) 



2/ t [l + Ml + ^ 2 )/2] 
j-(ZU Lp + NU Ln ) + 0^ F ) 



ZU Tp + NU Tn 
2«[l + ^(l + V 2 )/2] 

^-(ZU Tp + NU Tn )+0(C F ), 



TJ - — 



Gl p , n (r) + W 2p , n (r)A 



(4.14) 

(4.15) 
(4.16) 

(4.17) 
(4.18) 

(4.19) 



U Tp , n = 2tG 2 m ^(t) + W 2p , n (r)A , 



with 



and 



W 2p , n (r) 



1 + r 



G Ep,n( T ) + TG Mp,n( T ) 



T 
72 



i (e 2 F + e F yTTe + 1 + C 2 ) + A [e F + + + A 2 



(4.20) 



(4.21) 



-(1 + r). (4.22) 



Dividing the longitudinal and transverse RFG response functions by Gx and G^ yields the 
reduced responses 



^l RFG W = F^W = 



RFG, 
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2m N r] F 



;i-^(i-^ 2 



i + ^(i + ^ 2 ) 



(4.23) 



that, by construction, scale with ip, ( or y^ F G- 

In parallel with the scaling behavior of the RFG one can study the CSR and the vari- 
ous energy-weighted moments of another reduced response denoted r^K, A), introduced in 
Refs. [47,39]. Here the longitudinal response Rl(k, A) is divided by a function H l (k, A) to 
yield 



r L (ft, VO = \)/H l (k, A) 

and the n i/l moment of the longitudinal response of the nucleus is given by 



At 

» (n) = /" d\X n r L (K,X). 



(4.24) 



(4.25) 
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In particular, the n = moment, S(°), is the CSR. In the case of the RFG the dividing 
function is 

3 

H l (k, A) = (ZU Lp + NU Ln ) I (d^/d\) . (4.26) 

Thus, upon dividing the charge response of the RFG by Eq. (4.26) one obtains the 
following reduced longitudinal response 

^^-s^ 1 -™ 1 -^ (427) 

which, by construction, fulfills the CSR in the non-Pauli-blocked domain, as can easily be 
verified. 

The two dividing functions (related to the scaling) and H L (related to the CSR) are 
linked according to [39] 

G ^=(^)fx i + U^^)/2 H ^ X) (428) 
' 1 1 ^ //i.(k.A) • OiCy). (4.29) 



2 \tJ \ 1 + A # 

When experimental data are reduced using the above dividing functions, they are seen to 
yield a CSR at high-g and to scale when plotted versus any of the scaling variables previously 
introduced. 

Finally we explore the scaling and CSR properties of nuclear models other than the RFG, 
specifically: 

1) the hybrid model (HM), introduced in Ref. [39]; 

2) the quantum hadrodynamical (QHD) model [54]. 

The hybrid model is designed to account for the binding of the nucleons inside the 
nucleus, thus curing a flaw of the RFG related to its negative separation energy. The HM 
has continuum states that are plane waves, as in the RFG model, but has bound states 
described by shell-model wave functions obtained by solving the Schrodinger equation with 
some choice of potential well (in [39] harmonic oscillator bound-state wave functions were 
used to simplify the analysis in the limit where A — > oo). In the HM the scaling variable 
cannot be obtained analytically, since the model can only be treated numerically. However, 
scaling variables exist that incorporate the shift result from using Eq. (4.10) for (' (and 
a corresponding dimensionful variable y') or Eq. (4.11) for ip' by making the replacements 
A — > A' and r — > r' = k 2 — A' 2 , where 

A' = A - A shift , (4.30) 

with 

Ashift= 2^ {Tf + Es) (431) 
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FIG. 25. The reduced response F^ M (the hybrid model scaling function) is shown as a function 
of tp in panel (a) and ip' in panel (b) for 4 values of q (fine dotted — 0.5, dot-dashed — 1.0, heavy 
solid — 1.5 and solid — 2.0 GeV/c). The RFG result, which scales exactly as a function of tp, is 
also shown for reference as a dashed curve in panel (a). 

and T F = m N £ F the Fermi kinetic energy. The HM turns out to have the width of its 
reduced response identical to a RFG computed with a Fermi momentum that is somewhat 
larger than the usual one (237 MeV/c for the HM, versus the 230 MeV/c value for hp used 
for the RFG to correspond to nuclei near 40 Ca or 56 Fe). 
Scaling may then be examined for the HM by computing 

F L HM (/^) = R? m (k,\)/G l (k,\) (4.32) 

and the various energy-weighted moments of the longitudinal response, including the zeroth 
moment or CSR, may be computed using 

r™(K,i;) = R™(K,X)/H L (K,\). (4.33) 

Note that the same dividing factors Gl and that were developed from our discussions of 
the RFG are used. 

In Fig. 25 the scaling function F^ M for 40 Ca is displayed versus ip and ip' for four different 
values of q (the RFG result is also shown for reference): the HM scales either with ip or with 
ip' as q becomes large. Indeed, only the q = 500 MeV/c plot versus ip shows any appreciable 
violation of scaling, whereas the scaling versus ip' is excellent. 
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In the QHD model [54] protons and neutrons in the nucleus are described by Dirac 
spinors and move in strong Lorentz scalar and vector mean fields. These in turn arise self- 
consistently from the exchange of a and oo mesons between the same nucleons on which 
they act. The scalar field dresses the bare mass of the nucleon, considerably lowering its 
value; the vector field uniformly shifts the fermion spectrum. As a consequence the QHD 
charge response of nuclear matter in Hartree approximation is unaffected by the vector field, 
while it turns out to be quite sensitive to the effective mass m* N induced by the scalar field. 
This is, of course, true in the simple approximation of constant relativistic mean fields. An 
improved description allows for an energy- dependence of the latter, which helps to account 
for the data of proton-nucleus elastic scattering. 

As for the HM model, scaling can be examined in the QHD model by computing 

F L QHD ( K , V) = itf^, A)/C?l(«, A) (4.34) 

and likewise the various energy-weighted moments of the longitudinal response computed 
using 

r? HD ( K , V) = R? UD (k, X)/H l (k, A). (4.35) 

Since the dividing factors are (at least to a very good level of approximation) universal, 
accordingly we use the same dividing factors and H L that were developed from our 
discussions of the RFC 

In Fig. 26 we display the reduced responses in Eqs. (4.34) and (4.35) as functions both 
of ijj and also ip*, namely, the RFG scaling variable given in Eq. (4.11) with replaced 
by m* N (two different values of the effective mass, m* N = 0.68 m^v and 0.8 itln are used). It 
is clearly seen that _pQ HD does no t scale versus ip when the effective mass is constant and 
differs from m N . As q continues to grow beyond the range of values shown in the figures, 
the results continue to shift to higher to and never coalesce into a universal curve. When 
plotted versus ip* the behavior, while better, still does not scale. This is in contrast with the 
RFG and HM results displayed above and, importantly, is not what is seen experimentally 
where the world data do appear to scale in ip [51]. The fact that experimentally the scaling 
is observed to occur successfully for q greater than about 1 GeV/c suggests that m* N /m N 
should not deviate appreciably from unity for such kinematics. 

Finally, the CSR (S*- -*), the energy-weighted sum rule (S < - 1 - ) /S < -°- ) ) and the variance a = 
\JeS?> — (HW) 2 of 40 Ca corresponding to the HM, QHD (with m* N = 0.68 m N and 0.8 m N ) 
and RFG models are displayed in Fig. 27. The RFG model and the HM both saturate the 
CSR at high-g. In contrast, the QHD model does so only if the effective value of m* N /m N 
evolves with increasing q towards unity, as suggested by the latest version of the model. 



54 




FIG. 26. The reduced response with m* N = 0.68 mjy is shown as a function of ip in panel 

(a) and ijj* in panel (b) for 4 values of q (fine dotted — 0.5, dot-dashed — 1.0, heavy solid — 1.5 
and solid — 2.0 GeV/c). The RFG result, which scales exactly as a function of ip, is also shown 
for reference as a dashed curve in panel (a). 



V. OUTLOOK AND PERSPECTIVES 

In this paper we have discussed the parity-conserving and violating- inclusive nuclear 
responses, addressing the following issues 

a) Where do they occur? 

b) How can we describe them? 

c) What is their input? 

With respect to the first item we have limited our focus to the quasielastic peak (QEP), 
since at larger excitation energies where for example the A plays a role we do not expect 
a description of the nucleus only in terms of nucleonic and mesonic degrees of freedom to 
be adequate. At some point QCD degrees of freedom should become the more appropriate 
ones to describe inclusive scattering. On the other hand at energies significantly lower than 
those characterizing the QEP where discrete excitations, giant resonances, etc. are seen the 
theoretical many-body framework is different from the one used here. 

Concerning the second item let us again stress that any theoretical framework should 
first fulfill (as much as possible) Lorentz covariance and gauge invariance. For this the RFG 
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FIG. 27. The Coulomb sum rule (a), energy- weighted sum rule (b) and variance (c) are shown 
as functions of q for three models: dashed — the RFG, solid — the HM and the QHD model with 
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appears to be a good starting point: It is a covariant model, since its ingredients are the 
fully relativistic nucleon propagators and EM (or WNC) vertices, and it respects gauge 
invariance, since the vector currents of the nucleon are conserved. When mesons are added 
to the picture then Lorentz covariance and gauge invariance are fulfilled to the extent that 
one allows for a dynamical propagation of the mesons and treats the forces and the currents 
consistently. This turns out to be possible for the case of the pion. 

Of course the RFG misses surface and finite-size effects; however, first of all, these are of 
minor relevance for the scattering of electrons in the QEP and A peak domains, and secondly, 
they can be satisfactorily accounted for within the semiclassical approach, which exploits 
the advantages offered by the translational invariance of the RFG and accommodates these 
advantages to fit the physics of finite systems [56]. 

In the framework of the RFG we have treated nucleon-nucleon correlations in a pertur- 
bative scheme. Alternatives are represented by 

i) variational approaches 
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ii) loop expansions in the path integral framework, which however represent an alternative 
regrouping of perturbation theory. 

Concerning perturbation theory we have seen that basic ingredients are the diagrams 
associated with the HF mean field and those related to the antisymmetrized RPA. The 
emphasis on the antisymmetrization arises from the recognition that very important carriers 
of the nuclear force, notably n and p mesons, only act through the exchange diagrams. 

The short-range correlations (SRC) induced by the violent repulsion present in the N-N 
force at small distances are then inserted via the ladder diagrams. It is, however, doubtful 
whether ladders can be covariantly computed, especially at high densities (or Fermi mo- 
menta hp) where they become increasingly important and the role of relativity cannot be 
ignored. In addition, it may be that searching for a totally deterministic account of the N-N 
correlations is not the optimal way to proceed — should, for example, quantum chaos be 
at work in atomic nuclei, then costly efforts to compute SRC would not represent the most 
efficient way of interpreting the nuclear response functions. 

A further serious challenge facing the perturbative approach relates to the estimate of 
the size of diagrams that are not considered. In this connection it appears that only the loop 
expansion offers a consistent criterion (the number of loops) to organize the perturbative 
series in classes of homogeneous diagrams allowing at the same time an estimate of the 
rate of convergence of the resulting new expansion. However, although this is theoretically 
established, its practical implementation is far from trivial. 

Finally a few words on the input to be fed into the perturbative scheme are appropriate. 
In general the Bonn potential appears as a well-founded representation of the N-N force. 
Indeed it provides an excellent representation of the energy behavior of the N-N phase 
shifts in free space and of the deuteron's properties. In addition, it can be derived from 
a Lagrangian defined in terms of nucleonic and mesonic fields, and therefore cast in the 
framework of an effective field theory that, from the theoretical viewpoint, represents a highly 
desirable feature. In particular, it allows for a test of the gauge invariance of the theory, 
i.e. the consistency between the forces and currents. In this connection it is worthwhile to 
emphasize that the most direct way to assess the validity of the mesons-plus-nucleons model 
of the nucleus ultimately rests on a deeper understanding of the role of meson exchange 
currents (MEC) in the nuclear responses. Much has been done in this field, although much 
remains to be done. 

It seems to us that at present the analysis of the nuclear responses is best performed on 
the basis of a hadronic approach such as the one we have pursued using the Bonn potential. 
It is important to remember, however, that such potentials are being used for nucleons within 
the nucleus and that these are differently off-shell from the conditions found in studying N-N 
scattering or the ground state of the deuteron. 

In summary, in Section II of the present article we have analyzed the quasielastic response 
functions for inclusive electron scattering in the standard framework of HF plus RPA. In 
addition, the impact of short-range correlations on the effective particle-hole force has been 
explored within the context of the G-matrix approach. The analysis has then been extended 
in Section III to include the observables that occur in studying parity-violating electron 
scattering, with emphasis placed on the "new" nuclear axial response. In Section IV the 
attention has been directed towards several general, important properties of the response 
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functions, in particular their Coulomb sum rule and scaling/ superscaling behavior. Finally, 
a short account of what lies ahead has been provided in Section V. 



APPENDIX A: ELECTROMAGNETIC FORM FACTORS OF THE NUCLEON 

Here we give the formulae for the EM form factors introduced in the definition of the 
response functions in Eq. (2.8). 

1. Non-relativistic Fermi gas 

In a non-relativistic calculation one defines, in terms of Sach's form factors, 

A* = 2 (Ala) 

# )2 = 2rG£ )2 , 1 = 0,1, (Alb) 

where r = \Q 2 \/Am 2 N = (q 2 - uj 2 )/4m 2 N , G$ = G Xp + (-l) 7 Gx n (X = E,M). A typical 
parameterization (although many others are possible) is the dipolar-plus-Galster one, namely 



Ge p (t) = G v D (r) 
Gm p (t) = ti p G v D {r) 
" r) 

G En (r) = -^ n rG v D (r)Ur). 



G Mn ir) = » n G v D (T) (A2) 



Here G^(r) = (1 + A^r)~ 2 is the vector dipole form factor, with = 4.97, whereas 
fx p = 2.793 and ji n = —1.913 are the proton and neutron magnetic moments, respectively. 
For Ge„ we have adopted the Galster parameterization [55] with £„(t) = (1 + A n r) _1 and 



A n ^ 5.6. 



2. Relativistic Fermi gas 

In Ref. [32] it was shown that in an RFG calculation a very good approximation to the 
exact treatment of the EM vertices can be obtained with the following definitions: 



AD 2 _ 
Jl — 

Ai) 2 _ 
Jt — 



1 + r 



Gg } +tG 



2t 



1 + Ld/2rriN 



(i) 

M 

r {I) 2 
\2 M > 



2 k 2 



2m 2 



N 



(A3a) 
(A3b) 



where ip r is the relativistic scaling variable in Eq. (2.12), kp is the Fermi momentum and 
the other quantities have been defined in Eqs. (A2). 



58 



APPENDIX B: FIRST-ORDER SELF-ENERGY 



Here we give the analytic expressions for the first-order self-energy based on the potential 
in Eqs. (2.1)-(2.2). T,^(k) is the sum of direct ("Hartree") and exchange ("Fock") terms, 
namely 



where 



and 



S (1) (A;) = £ H (A;) + £ F (A;), 
S H (A;) = pV o (0) 



8 



(Bl) 
(B2) 

(B3) 



p = 2kp/3ir 2 being the nuclear density. In the last equation we have introduced the spin- 
isospin coefficients (note that the tensor channels do not contribute) 



4 0) 



i, 4 T) = 

cf } = 0, 



3, a 



if) 

F 



3, Ct ] =9, 



and defined 



(B4) 



(B5) 



In any non-tensor channel a the potential is expressed as a combination of the "<5" and 
"momentum-dependent" pieces in Eq. (2.2), for which one finds 



„ 2 
96v 



and 



'md(^) — * 



£ = 

Sj(k) = { g s (X 2 - p 2 )w F a (X\k), t=\ 
g s (\ 2 -p 2 ) 2 w F (\\k), t = 2 



' gMDV 2 wl(\\k), 
g M B^ 2 [w F (p\k) -wl(X\k)\, 



(B6) 



e = o 
e = i 



{ 9M^>Mk) - w F a (X\k) - (A 2 - p 2 H(\\k)}, i = 2, 



(B7) 



where £ represents the power of the form factors (see Eq. (2.2)). Here we have introduced 
the dimensionless form factor cut-off, A = A/k F , and meson mass, p = m/kp, and we have 
defined 



w F a (\\k) 



w F b (\\k) = - 



A 



arctan 



k " 



A 



+ arctan 



1 + k 



arctan 



1-k' 



+ arctan 



A 



A 2 - k 2 + 1 



In 



A 2 + (k - 1) 



4k 



4k 

\ 2 + (k-l) 2 



\ 2 + (k + l) 



A 2 + (k + 1) 2 



(B8a) 
(B8b) 
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APPENDIX C: TENSOR INTERACTION IN THE EXCHANGE DIAGRAMS 



The n-th order exchange polarization propagator in presence of tensor interactions has 
an expression that is slightly more complicated than that in Eq. (2.60), because the tensor 
operators do not in general allow for a factorization of the azimuthal integrations. A generic 
diagram with m non-tensor and n — m tensor interaction lines can instead be written as 



n+l / , \ 2n+2 



(\ n+L / t, \ ■ 



x / dyi\ I 1 dxi ■ ■ ■ [ dy n+ i\ f " +1 dx n+1 

J-1 Z JO J-i Z Jo 

x- : — W ai (xi,yi;x2,y2) • • • W am ( 

ip-yi + 

x W am+1 ... an (x m+ i, y m +i, ...; x n+ i, y n +i 



*/> ~ Vn+i + iriu 

+ J>^-a;), (CI) 
where W ai has been defined for the non-tensor channels in Eq. (2.59) and 



Wa m +i...a. 



I \ _ nn-m V- V- a (** d <fm+l f 2n #n+l 
n {X m +l, Um+l', Xn+1, Vn+l) — z 2.^ 2.^ / — o "■ / 

„■■ , , Jo Ztt Jo 



2ir 



ij h...l„ 

x V am+1 (k m+ i — fc m+2 )S'j; 1 (k m+ i — k m+ 2)---V an (k n — k n+ i)Si n _ m j{k n — k n+ \). (C2) 
In the last expression we have introduced the tensors 

Sij(k) = 3kikj — 5ij, (C3) 

such that <TiCTjSij(k) = S 12(h). 

The first-order case is rather simple, since one again gets Eqs. (2.61)-(2.63) with 

W a (x,y;x',y') = [*" ^V a (k - k')S zz (k^k>). (C4) 
Jo Zn 

At second order, however, one can use Eqs. (2.64)-(2.66) only when just one tensor interac- 
tion is present. 



APPENDIX D: FIRST- AND SECOND-ORDER EXCHANGE DIAGRAMS 



Here we give the explicit expressions for the first- and second-order exchange diagrams, 
based on the potential in Eqs. (2.1)-(2.2). In Eqs. (2.61) and (2.62) we have seen that 



U^(q,cv) 



m N 



2 VA 



(2tt)< 



Qi 1} (0, 1>) - Q^(q^) + Q^M + q) - Qg\-q, 1> + q) 

(Dl) 
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where 



Qg\q,f/>) = 2 /* dy- 1 —— f dy'W a "(y,y';q) L— , (D2) 

J-i ip-y + irju, J -i y-y' + q 

whereas from Eqs. (2.64) and (2.65) one has 



q J (2tt) 6 

Qi 2 i(0, 0;^ + q) + 2^(0, -g; V + g) + Q&(-?, 0;^ + q) - G&(-?, "95 V> + , 

(D3) 



where 



(2) _ Z* 1 1 f 1_y2 _1 

Qaa>(QuQ2\^)= dy- dxG a (x,y + q^ip + qi)- — — Q a ,{x, y + q 2 ; ip + q 2 ) 

J-i 2 Jo ib — y + 



ip-y + irj u 

(D4) 



and 

- 1 . , 1 



Q a {x,y;i>)= I dy'- ——W^x,y;y'). (D5) 

y-i ip-y' + n}u 

For a meson-exchange potential the quantities that can be calculated analytically are those 
given by Eqs. (2.59), (C4), (2.66b) and (2.63), namely 



W a (x,y;x',y') = f -^V a (k - k') (non-tensor) (D6a) 
Jo 2ir 

W a (x,y;x',y') = [*" ^V a (k - k')S zz {k^k') (tensor) (D6b) 
jo 2ir 



and 



W'Jx, y; y') = \ f * dx' W a (x, y; x', y') (D7) 
2 Jo 

W^(y, y'; q) = \ dx \ f~ V dx' W a (x, y + q- x', y'). (D8) 

2 Jo 2 Jo 



In any channel a the potential is expressed as a combination of the terms displayed in 
Eq. (2.2). Then, for each of them one finds 



W s (x,y;x',y') = < 



96, i = 

g s (X 2 - fi 2 )w a (X\x,y;x',y'), £=l (D9a) 

[ g s (X 2 -^ 2 ) 2 w b (X\x,y;x',y'), £ = 2 

gMB^ 2 w a (^\x,y]x',y'), £ = 

fr ,,.„.> j. ;^MD^ 2 K(/x|x,y;x',y') - w a {\\x,y;x' ,y% £=l , . 

W MD (x, y,x,y)-{ gMD AwM ^ y _ x , yl) _ y _ x , yl) (D9b) 

-(X 2 -^ 2 )w b (X\x,y;x',y% £ = 2 
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W TN (x,y;x',y') 



' ^TN{[3(y - y'f + fx 2 ]w a (fi\x, y; x', y') - 1}, £ = 

0tn{[3(2/ - y') 2 + n 2 ]w a (n\x, y; x', y') 

-[3(y-y'f + X 2 }w a (X\x,y;x',y')}, € = 1 

^TN{[3(y-y') 2 + ^] 

x [w a (n\x, y; x\ y') - w a {X\x, y; x\ y')\ 

-(A 2 - ^) [3(y - y'f + A> fe (A|x, y; x', y')}, £ = 2, 



(D9c) 



where again £ labels the power of the form factors, we have introduced the dimensionless 
form factor cut-off, A = A//cp, and meson mass, /x = m/kF, and we have defined 



w a (X\x, y; x', y') = {[A 2 + (y - y'f + x + x'f - Axx'}^ 2 

A 2 + (y - y'f + x + x' 



w b (X\x,y;x ,y ) 



{[A 2 + (y- y'f + x + x'} 2 - 4xx'} 3 / 2 ' 



(DlOa) 
(DlOb) 



For W' one finds 



W' s {x,y-y') 



WM D (x,y;y') = 



W^(x,y;y') = < 



g 5 (l-y' 2 )/2, £ = 

g s (X 2 - n 2 )w' a (X\x,y;y'), £=1 
9s(X 2 - fi 2 fw' b (X\x, y; y'), £ = 2 

' g M D^ 2 w' a (fi\x,y,y'), 
9md /i 2 K(/i|x, y; y') - w' a (X\x, y; y% 
g M DV [w' a (fi\x,y;y') - w' a (X\x, y; y') 
-(A 2 -^ 2 )w' b (X\x, y;y% 

[ <? TN {[3(y - y'f + vVMx, v; y') - (i - y' 2 )/2}, 
^TN{[3(y - y') 2 + A w M x i y; y') 

-[3(y-y'f + X 2 ]w' a (X\x,y;y')}, 
£tn{[3(?/ -y'f + fi 2 ][wMx,y;y') - w' a (X\x, y; y')\ 
-(X 2 -/, 2 )[3(y-y'f + X 2 ]w' b (X\x,y;y')}, 



£ = 
1 = 1 

£ = 2 



(Dlla) 



(Dllb) 



£ = 
£= 1 
* = 2, 



(Dllc) 



where 



™l(A|x,y;y') = - 



x In 



w' b (X\x,y;y') = \ 



A 2 + (y - y') 2 + 1 - y' - x + J[X 2 + (y - y') 2 + 1 - y /2 + a;] 2 - 4(1 - y' 2 )x 



(D12a) 



4A 2 + (y- 2/ ') 2 



2[A 2 + (y-y') 2 ] 

A 2 + (y-yQ 2 -l + y /2 + x 

^[A 2 + (i - y'f + l-y' 2 + xf - 4(1 - y' 2 )x_ 

(D12b) 



Finally, for W 7 "" one finds 



r^[(l-y 2 )/2][(l-y' 2 )/2], £ = 
= y 5 (A 2 -/i 2 K'(A|y,y';g), €= 1 
(g s (X 2 -^fw' b \X\y,y';q), £ = 2 



(D13a) 
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g M D^ 2 wMy,y';q), £ = 

WM D (y,y';q) = { g M D fi 2 [wMy,y';q) - w'^(X\y,y';q)], £ = l 

g M D^>My,y';q)-<(My,y';q) - (x 2 - M(\\y, y f -q)}, £ = 2 

(D13b) 

f 9TN{[3(y -y' + qf + fi 2 }wMy,y>; q) 

-[(l-y 2 )/2][(l-y' 2 )/2]}, £ = 

9TN{[3(y ~ y' + q) 2 + /i 2 K'(/i|y,y';g) 

-[3(y ~y' + qf + A 2 K(A|y, y'; q)}, £ = 1 imdCj 
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^ N (y,y';?) 
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